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Der Beitrag von weicher Streuung zum Jet Quenching in einem Quark-Gluon 
Plasma und allgemeine Eigenschaften des Energieverlusts von Partonen 

Ausgehend von der Beobachtung, dass weiche Streuungen eines hart gestreuten 
Partons am Farbladung tragenden Protoniiberrest in Proton-Proton Kollisionen 
eine wichtige Rolle spielen, wurde ein ahnliches Modell fiir weiche Streuung in 
einem Quark-Gluon Plasma konstruiert und als Monte Carlo Generator implemen- 
tiert. Die hauptsachliche Frage war, inwieweit diese weichen Streuungen zum En- 
ergieverlust beitragen konnen, den ein hochenergetisches Parton erleidet, wenn es 
ein Quark-Gluon Plasma durchquert. Das Ergebnis war, dass sie bis zu 50% des 
beobachteten Energieverlusts erklaren konnen. 

AuBerdem wurde die Frage untersucht, was man von heute verfiigbaren Daten 
liber die allgemeinen Eigenschaften des Energieverlusts von Partonen erfahren kann. 
Dazu wurde ein Monte Carlo Modell benutzt, das zwar die voile Simulation des 
Quark-Gluon Plasmas beinhaltet aber fiir den Energieverlust einen allgemeineren 
Ansatz verwendet. Dadurch wird es moglich verschiedene Szenarien wie koharente 
Gluonbremsstrahlung, Streuung etc. nachzubilden. Leider stellte sich heraus, dass 
es sehr schwierig ist die verschiedenen Moglichkeiten zu unterscheiden, sodass zum 
jetzigen Zeitpunkt keine eindeutige Entscheidung zugunsten der einen oder anderen 
moglich ist. 



The Soft Scattering Contribution to Jet Quenching in a Quark-Gluon 
Plasma and General Properties of Partonic Energy Loss 

Starting from the observation that soft rescattering of hard scattered partons from 
the colour charged remnant plays an important role in proton-proton collisions a 
similar model for soft scatterings in a quark-gluon plasma was constructed and 
implemented as a Monte Carlo event generator. The main emphasis was put on 
the question to what extend these soft scatterings can contribute to the energy loss 
that an energetic parton suffers, when it traverses a quark-gluon plasma. It was 
found that the soft scattering can account for up to 50% of the observed energy 
loss. 

Furthermore, it was investigated what information present data reveal about the 
general features of partonic energy loss. This study was carried out with a Monte 
Carlo model that includes the full simulation of the quark-gluon plasma but uses 
a more general ansatz for the energy loss which makes it possible to emulate dif- 
ferent scenarios like coherent gluon bremsstrahlung, scattering etc. Unfortunately, 
it turned out that it is difficult to differentiate between them so that at present no 
clear decision in favour of one or the other is possible. 
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1 Introduction 



The quark-gluon plasma (QGP) is a state of deconfined quarks and gluons and re- 
stored chiral symmetry predicted by Quantum Chromodynamics (QCD), the theory 
of the strong interactions. It is expected to be formed at very high energy densities. 
It presumably existed in the early universe and can possibly still today be found in 
some kinds of neutron stars or other exotic objects. Great efforts are being made 
to produce the QGP in the laboratory by means of heavy ion collisions. This is the 
only possibility to reach high enough energy densities in man-made experiments. 
The price that has to be paid is that the QGP - if it is formed at all - has a very 
short lifetime because of rapid expansion and thus cooling of the system. 

The detection and study of the QGP is a highly non-trivial task that has kept 
both experimentalists and theorists busy for many years and is still far from com- 
pleted. The observables can be grouped into soft and hard probes. The hard probes 
are sensitive to the very early stage of the collision whereas the soft ones probe the 
latest stage. Hard observables such as jet quenching (suppression of particles with 
high transverse momentum) have only recently become accessible with the Rela- 
tivistic Heavy Ion Collider (RHIC) in Brookhaven, where data taking started in 
2000. They are considered an adequate tool for investigating the properties of the 
medium created in the collisions wherefore sometimes the term "jet tomography" 
is used. For the proof of QGP formation on the other hand other observables like 
thermal and chemical equilibration are more suited. 

Most of the phenomena connected to heavy ion collisions and QGP formation 
and its evolution belong to the non-perturbative regime of QCD where it is very 
difficult to derive a quantitative description from first principles. Thus models of 
different scenarios have to be constructed and the results compared to data in order 
to learn something about the physics. This is often a laborious business and it is 
in the majority of cases not possible to obtain clear cut conclusions. 

In this thesis a new model for the jet quenching is presented and compared to 
data. It is built on partonic energy loss due to soft colour interactions. The starting 
point was the Soft Colour Interaction (SCI) model [l] that successfully describes a 
wealth of data mainly on deep inelastic scattering and diffraction. The SCI model 
has been extended in order to be applicable for the jet quenching phenomena. 

There is no perfect model for the jet quenching, so this new model can hopefully 
help to gain a better understanding of the effects in QGP and soft physics in 
general. The SCI jet quenching model is implemented as a Monte Carlo event 
generator which has the advantage that the processes can be simulated and studied 
in great detail. 
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A second issue of this thesis are the general properties of the energy loss mecha- 
nism. Different theoretical approaches lead to quite different behaviours. The idea 
is to see what presently available data reveal about the main features of the sought 
energy loss mechanism without having a detailed model. This information may 
then help to get to the right description in the end. 

In the first part of this report a short overview of QCD, the search for the QGP, 
jet quenching as a possible signature of the QGP and a few basics of heavy ion 
physics is given. There is also a discussion of experimental results with the main 
focus on jet quenching. After a brief introduction to Monte Carlo methods and the 
SCI model the new jet quenching model is described and compared to data in more 
detail. In the following section an investigation of general features of energy loss in 
the QGP is carried out with the help of an adjusted Monte Carlo model. 
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2 QCD for Beginners 



This section gives an introduction to the basics of QCD and some topics that are 
relevant for this study. A more systematic and detailed discussion can be found in 
12 0] or other standard books. 

2.1 Introduction 

Quantum Chromodynamics (QCD) is the theory of the strong interaction. It de- 
scribes the interaction of quarks via exchange of gluons. There are the six quarks 
called down (d), up (m), strange (s), charm (c), bottom (b) and top (t) which are 
fermions and the massless gluons which are bosons. The relevant charge, which is 
the source of the field, is called colour and the three states are labelled red (r), blue 
(6) and green (g). Quarks come as colour triplets, which means that a particular 
quark can carry any of the three colours. The antiquarks form the corresponding 
antitriplet (i.e. they carry the anticolours f, b and g). Quarks cannot be observed 
as free particles but are always confined into hadrons, which are colour neutral ob- 
jects. A colour neutral combination is also called a singlet. There are at least two 
possible neutral colour combinations: a colour and its matching anticolour (giving 
rise to mesons that are composed of one quark and one antiquark), or the three 
(anti)colours (leading to (anti)baryons that consist of three (anti) quarks). There is 
also a possibility to combine two quarks, but it is not clear to what extend these 
diquarks are bound states. Nevertheless, they can be very useful, e.g. for the treat- 
ment of remnants as will be discussed later. The combination of two triplet charges 
gives an antitriplet (or a sextet, but this is not interesting in this context), two 
antitriplets form a triplet. 

The gluons are also colour charged, but unlike the quarks they carry octet charges 
which can be viewed as the combination of a colour and an anticolour (the com- 
bination of colour and anticolour gives a singlet and an octet). The name octet 
already indicates that it contains eight states. The naively expected nineth gluon 
does not occur because as a colour singlet combination it does not correspond to 
an interaction. The fact that the gluons carry colour charge gives rise to the gluon 
self-interaction, which means that gluons can directly couple to each other. 

The coupling strength of the strong interaction is not constant but depends 
on the momentum transfer such that it increases with decreasing energy scale. 
Therefore processes involving large momentum transfers can be treated analytically 
since it is possible to make a perturbative expansion in powers of as- In regimes with 
small energy scales, where the coupling is large, pertubation theory is not applicable 
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and other ways to treat these problems have to be found. One possibihty is to 
solve the equations numerically, which has its own problems. This task is highly 
non-trivial and to a great extent limited by the available computing power. In 
many cases one has to rely on phenomenological models that describe the observed 
phenomena without having a firm theoretical basis. 



2.2 Structure of the proton and parton showers 

The proton is a baryon built of three valence quarks (two u- and one d-quark) 
which carry its electric charge and baryon number. In addition, there are also 
gluons and sea quarks that are fluctuations of gluons into quark-antiquark pairs. 
The quarks carry approximately 50% of the proton's energy-momentum, the other 
half is carried by the gluons. 

fc^^^^^^^^ k' 



Figure 2.1: Inelastic lepton-proton scattering via exchange of a virtual photon 

The best way to investigate the structure of the proton is in lepton-proton scat- 
tering via photon exchange fFig. 12. ij) . The standard variables are 



Q' = -q' (2.1) 

X = — ^ (2.2) 
2p-q ^ ^ 

y = k-p ^2.3) 

where p, q and k are the proton's, the photon's and the lepton's four-momenta, 
respectively {k refers to the initial and k' to the final vector). 

At small momentum transfers the photon interacts with the proton as a whole, 
but at higher the structure of the proton can be resolved. The momentum 
transfer can be viewed as the resolution: With a bigger smaller structures can 
be resolved (the resolution is (i ~ 0.2fm/ y^Q^[GeY]). The cross section can be 
written as 

-iil-y)F, + xy^F,) (2.4) 



dQ^dx xQ^ 
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where the structure functions Fi^x^Q"^) parametrise the proton structure. They are 
approximately independent of (Bjorken scahng) which imphes that the photon 
interacts with quasi-free pointhke constituents (partons) of the proton. In the frame 
where the proton has infinite momentum, x can be interpreted as the hmit of the 
fraction of the proton's momentum carried by the parton. The structure function 
F2 can to lowest order be written as 

F2(x,g2) = ^e2x/,(x,g2) (2.5) 

a 

where a runs over all flavours and antiflavours, is the electric charge of the 
respective quark or antiquark and fai^iQ'^) is the probability of finding a quark of 
flavour a carrying the momentum fraction x inside the proton when probing with 
momentum transfer Q^. The fa are called parton distribution functions (pdf's), 
they include the contributions from both the valence and the sea quarks. Scattering 
of the lepton from gluons is not possible since the gluons do not carry electric charge. 
Fi is connected to F2 via the Callan-Gross relation 

F2 = 2xFi (2.6) 

The dependence of the parton distribution arises from the fact that the 
(anti) quark may have radiated gluons (initial state radiation) before it was ac- 
tually struck by the photon. The dependence can be calculated in pertubation 
theory and is governed by the DGLAP^ evolution equations 



dlnQ^ 271 



dy 

y 



X 



19 



X 



fg{y,Q' 



(2.7) 



dfg{x,Q') _ as{Q' 



dy 

y 



P, 



91 



X 



y 



f,{y,Q' 



P. 



99 



y 



(2. 



where fq and fg are the quark and gluon pdf's, respectively. The Pij are called 
splitting functions and describe the evolution of a parton j into a parton i carrying 
the energy fraction z = x/y oi the original parton. Pqq for instance describes the 
radiation of a gluon from a quark where the quark keeps the energy fraction z\ 
Pgq describes the same situation when the gluon gets the fraction z and the quark 
obtaines 1 — z. 

There are similar equations for gluon radiation after the photon-parton scattering 
(final state radiation). The gluons tend to be emitted with small angles relative to 
the (anti) quark. The branching processes can explicitly be simulated using parton 
showers. 



^Dokshitzcr Gribov Lipatov AltarcUi Parisi 
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In contrast to the Q'^ dependence, the x dependence of the pdf's has to be 
parametrised and fitted to data. There are of course parton densities for all hadrons, 
it is common to write fajh for the pdf of the parton a in hadron h. 

For (perturbatively hard) hadron-hadron processes any inclusive cross section 
can be factorised into the form 




where (Ja^h is the parton level cross section for the parton species a, 6 and ^ is 
an arbitrary factorisation scale separating the hard from the soft regime. (It is a 
factorisation of the process in pdf's and the parton-parton process described by 



2.3 Fragmentation 

Quarks and gluons cannot be observed as free particles but are confined in colour 
neutral hadrons. The process of hadron formation from partons (quarks and gluons) 
falls into the domain of non-perturbative QCD and it has not been possible to derive 
it from first principles (it might become possible with lattice QCD). Nevertheless, 
a profound understanding of these mechanisms is of great importance for all kinds 
of activities in high energy physics, where the dynamics on the parton level has to 
be deduced from the observed hadronic final state. 

There are basically three models: string fragmentation, independent fragmenta- 
tion and cluster fragmentation. The main ideas of the Lund string model j3] and 
the independent fragmentation ([S], see also jH] for an overview) will be reviewed 
briefly here. 



2.3.1 String Fragmentation 




Figure 2.2: Quark-quark scattering via gluon exchange as an example of inelastic 
pp scattering with the resulting colour field topology represented by 
strings (dotted lines) 
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The jet production from a hard quark-quark scattering in a proton-proton col- 
hsion can serve as a simple example for the fragmentation of a quark. In this 
interaction the protons are broken apart: The two quarks are scattered to large 
angles and the remnants (in this case diquarks) continue in the beam directions. 
The incoming protons are colour singlets, after the scattering the (colour-) charged 
constituents can again be grouped into two colour neutral systems (Fig. I2.2|l . A 
field stretches between a colour charge (triplet charge) and its matching anticolour 
(antitriplet, in this case the diquark). Due to the gluon self-interaction the field 
does not extend transversly in space but can be viewed as a colour flux tube with 
transverse size ~ 1 fm. The potential rises linearly with the distance between the 
charges and at some point it becomes possible to form a new quark-antiquark pair 
from the field energy. This process is sketched in Figure 12. HI 




q qq q qq q q q qq 



(a) (b) (c) 

Figure 2.3: Creation of a qq pair from the field energy during the separation of a 
quark (g) and a diquark {qq) 



In the Lund model the colour flux tube is idealised to a one-dimensional massless 
relativistic string with energy density k ~ IGeV/fm. A newly produced qq pair 
breaks the string into two independent parts. The creation of a qq pair is described 
as a quantum mechanical tunneling process so that the probability becomes jlj 




(2.10) 



This leads to a Gaussian distribution of the two transverse components of the 
produced (anti) quark's momentum {p^^ and Py if the string is stretched along the 
z-axis). The transverse momentum is compensated when the string breaks between 
the created quark and the antiquark. Equation 12. 101 shows also that the creation of 
heavy flavours is suppressed so that charm, bottom and top quarks are not expected 
to be produced. Since the quark masses are not well defined in this context the 
suppression of strangeness relative to u (or d) quarks is left as a parameter with 
7, = ^ ~ 0.3. 

It is assumed that a meson is formed when the invariant mass of a string piece 
becomes small. Its p±_ is given by the sum of the transverse momenta of the string 
endpoints, which also fix the flavour composition. The spin of the meson is chosen to 
be either or 1 according to a suitable probability, orbital excitations are expected 
to be rare and are therefore excluded in the default Lund model. 

The different break-ups of the string are assumed to be independent of each other. 
It is therefore possible to start at one end of the string, make a break-up and be 
left with a meson and a shorter string. The procedure can be iterated until the 
available energy is used up (the termination requires some extra treatment that will 
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not be discussed here). Let us consider the string in the cm. frame of the quark 
and the diquark with the quark moving in the +z-direction and start with the 
fragmentation from the quark end. The first step is to choose the flavour, spin and 
transverse momentum of the created qq pair as described above. The identity and 
transverse mass of the meson are then fixed and only the longitudinal momentum 
(or the energy) remains to be determined. This is done by assigning the meson a 
fraction z of the available E + p^. What is left to the new string is j3] 

(E + p,)ncw = {l-z){E + p,Ud (2.11) 

9 

(E-p,)new = (^-p.)old- , ^ X (2.12) 

The values of z are distributed according to a probability distribution f{z). The 
constraint that the result should be independent of the choice from which end to 
start leads to the 'Lund symmetric fragmentation function' 

f{z)ocz-\l-zyexp(^-^^ (2.13) 

One important point is that the hadron that contains the end-(di)quark of the 
string will typically have the highest momentum in the ensemble. It should further 
be noted that in the lab system it is the slowest hadrons that are formed first, 
because all particles have a formation time of ~ 1 fm in their rest frame. The 
difference in formation time in the lab frame is a boost effect. This ordering is 
frame dependent and therefore not important for the fragmentation algorithm, 
which is Lorentz invariant. 



qq , 




Figure 2.4: Inelastic pp scattering with additional gluon radiation resulting ac- 
cording to the Lund model in a colour string field topology (dashed 
lines) where radiated gluons are connected in the same string as the 
scattered quarks 



The formation of baryons is somewhat more complicated but works in principle 
like the meson production. In the 'popcorn model' baryons arise from cases where 
the produced qq pair does not match the colour of the string ends. It is then possible 
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that another pair with the third colour is produced and a baryon and an antibaryon 
is formed. Baryon production is suppressed since pairs with the "wrong" charge 
can only exist as fluctuations and also because two qq pairs must be produced. 
Another source of baryon production in the Lund model are the remnant diquarks, 
which are turned into a baryons via normal string break-ups with formation of a 
qq pair. 

Gluons carry colour octet charges which means that they cannot serve as endpoint 
for a triplet string. But they can be situated in the middle of a string with the 
colour connected to the anticolour-end and the anticolour with the colour-end of 
the string. If a quark for example radiates gluons in a parton shower a la DGLAP 
the gluons will be aligned in the same string as the quark (Fig. I2.4|l . 

2.3.2 Independent Fragmentation 

In the framework of independent fragmentation it is assumed that each parton 
hadronises on its own, i.e. the fragmentation of a jet system is an incoherent sum 
of the fragmentations of each parton. The procedure has to be carried out in the 
overall cm. system. 

There is an iterative approach somewhat similar to the Lund model: A jet arising 
from the fragmentation of a quark is split into a meson and a remainder jet with 
lower energy and momentum. The sharing of energy-momentum is governed by 
a fragmentation function, actually the same functions can be used for string and 
independent fragmentation. Flavour and transverse momentum are conserved in 
each break-up. It is assumed that the formation of a meson does not depend on the 
energy of the remainder jet so that the step can be iterated resulting in a sequence 
of hadrons. There is, however, one problem with very small values of z: They lead 
to backward moving hadrons (i.e. < 0) that have to be rejected. 

Although flavour and transverse momentum are conserved locally this is not the 
case in the global balance. In the end there will always be an unpaired (anti) quark 
left behind. There may also be hadrons with p\\ < 0, that were removed their 
energy and momentum being lost for the jet. Thus overall energy, momentum, 
charge and flavour are not conserved. This also applies to jet systems since each 
parton hadronises separately. 

There are several possibilities for the treatment of gluon jets. Since a gluon is 
expected to result in a softer jet it is sensible to split the gluon perturbatively in a 
qq pair, then the same fragmentation function as for the quark jets can be used. 

Apart from the non-conservation of energy-momentum and flavour there are two 
more conceptual weaknesses: the issues of Lorentz invariance and coUinear diver- 
gences. The result of independent fragmentation depends on the coordinate frame 
and is thus not Lorentz invariant. This problem is circumvented by requiring the 
fragmentation to take place in the overall cm. frame. The colliniear divergence, on 
the other hand, leads to problems in connection with parton showers. A system of 
collinear partons leads to a much higher hadron multiplicity than a single parton 
with the same energy. 
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3 Hunting the Quark Gluon Plasma 



3.1 QCD Predictions 

At high temperatures and densities the long range interactions between quarks 
are dynamically screened (similar to Debye screening). Only the very short range 
interactions remain but here the coupling is weak so that the quarks and gluons 
are quasi-free and thus deconfined. Furthermore chiral symmetry is restored at 
apparantly the same critical temperature Tc. This phase of deconfined quarks 
and gluons and restored chiral symmetry is called the quark-gluon plasma (QGP) 
[7]. This is not unexpected since spontaneously broken symmetries (such as chiral 
symmetry) are often restored at high temperatures through a phase transition jS] . 
An important question is now whether there is a phase transition from hadronic to 
deconfined matter. Recent results of simulations of QCD on the lattice indicate a 
phase transition at a critical temperature Tc ~ 170 MeV which corresponds to an 
energy density of ec — 1 GeV fm~^ jHj . This is too low for pertubation theory to be 
applicable and one has to rely on lattice QCD. Although great progress has been 
made in this field there are still major problems. It has, for instance, so far not been 
possible to determine the order of the phase transition and many calculations are 
done for vanishing baryon chemical potential (i.e. vanishing baryon number) which 
is a good approximation for the early universe and the LHC but not for AGS, SPS 
and RHIC [71 E]. 

3.2 Models 

There is a large variety of different models for heavy-ion collisions that are based 
on largely different ideas and assumptions. Only a short overview over the main 
classes can be given here [3 |Hj • 

Statistical models (e.g. ^Oj) assume that local thermal and chemical equilibrium 
is achieved during the collision. The starting point is a hadron gas that can 
either be created directly in the collision or the product of the hadronisation 
of a QGP. It is described as an ideal hadron gas using a canonical or grand 
canonical formalism. The gas expands until inelastic interactions cease, then 
the composition of the system is fixed (chemical freeze-out). At some point 
the mean free path becomes so long that also the elastic interactions stop 
(thermal freeze-out). When resonance decays are included statistical models 
can yield the relative abundancies of hadron species as they are measured by 
experiments. 
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Comparison of the model results with data can help to clarify if, or to what 
extent, equilibrium is achieved. The freeze-out temperature and (hadron) 
chemical potential can be determined by fitting the model to data. 

Parton Cascades (e.g. [TT]) are detailed microscopic models that describe the col- 
lision of two nuclei in a perturbative QCD framework. The first step is the 
decomposition of the nuclei in partons according to measured structure func- 
tions. The interactions during the collisions are treated as perturbative scat- 
terings with initial and final state radiation. The last stage is the hadronisa- 
tion of the partons using the Lund string model. 

Parton cascades predict a rapid thermalisation (proper time scale 0.3-0.5 fm) 
and a chemical equilibration that takes somewhat longer (several fm). The 
plasma is initially gluon rich due to the larger cross sections for gluons. 

Hadronic Transport models (e.g. [T2j) are formulated in a hadron basis although 
some also include non-hadronic elements such as quarks and strings. The 
heavy-ion collision is described as a sequence of collisions of constituents 
(mesons, baryons, quarks,...). Partonic degrees of freedom are not treated 
explicitly so that no phase transition can occur. 

Hadronic models provide a very useful background for other models because 
it is important to understand which phenomena can be described in terms of 
hadronic physics. 

Hydrodynamic models (e.g. O Ej) are macroscopic kinetic models that are 
based on the assumption of local equilibrium and energy-momentum con- 
servation. The nuclei are described as fluids and in some models a third fluid 
can be created in the collisions. Starting from the colliding nuclei, an equi- 
librated QGP or hadronic matter the time evolution of the system can be 
studied until hadronic freeze-out. 

These are the only dynamical models in which a phase transition can be 
incorporated explicitly via the equations of state. 

3.3 Observables 

The number of observables that have been suggested as signatures of QGP forma- 
tion is so large that it is impossible to discuss all of them here. Instead three of the 
most popular will be presented without going into details. 

Strangeness Enhancement 

The production of strange hadrons is suppressed in pp colhsions and the suppres- 
sion increases with the strangeness content of the respective hadron. This has been 
argued to be due to the higher strange quark mass [3 IHj- In a QGP strangeness 
saturation via ss pair production is expected which would significantly increase 
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the yields of strange particles jH]. The estimated time scale of several fm for 
strangeness equilibration is maybe too long for a complete saturation, but an in- 
creased strangeness content leads to increased strange particle yields in any case. 
If it is the statistical hadronisation of a strangeness-enhanced deconfined phase 
that is observed a stronger enhancement of multistrange hadrons is expected. The 
enhancement factor for a hadron containing strange quarks is where Eg is 
the global enhancement factor. Hadron rescattering scenarios lead to the opposite 
behaviour p!3] . 

A QGP formed at AGS or SPS would have nonzero chemical potentials for u and 
d leading to an ordering in the quark abundances: The densities of u and d is higher 
than that of s and s, which is higher than the u and d density. This means that at 
freeze-out the combination of a s with a n or d to K'^ or is more likely than the 
su and sd combination {K~ and K'^). Thus the /ti^ ratio should be different 
from K~ /t:~ in case of QGP formation. Unfortunately, the argumentation has 
several drawbacks. One is that the strange particle abundancies after freeze-out of 
a QGP are very close to those in an equilibrated hadron gas with the same entropy 
content jZj so that it is difficult to unambiguously relate different ratios to QGP 
formation. 

Chemical Equilibrium and Freeze-Out 

If the observed hadrons are produced in the hadronisation of an equilibrated QGP 
they should inherit the property of chemical equilibrium. This also includes the 
disappearance of the strangeness suppression. From a statistical model fit to data 
the freeze-out temperature and chemical potential can be found. These values can 
then be compared to lattice QCD calculations for the phase boundary. A chemical 
freeze-out close to or at the phase boundary suggests that the hadrons originate 
from a deconfined medium and that the chemical composition is established during 
the phase transition. 

Charmonium Suppression 

In a QGP colour screening reduces the range of the attractive force between quarks 
and antiquarks and thereby prevents cc pairs from binding. It is therefore expected 
that the more loosely bound iIj' and Xc states start to be suppressed at a lower 
temperature than J /ip jTHj. There is evidence from lattice calculations that the IS* 
states J/ip and rjc survive up to 1.5 Tc whereas the Xc states are dissolved already 
at 1.1 [in]. 

cc pairs are produced as small configurations and the evolution to the larger 
charmonium state takes approximately 1 fm. Thus the J / ip should survive if it 
escapes fast enough, i.e. if it has a high transverse momentum or if the QGP 
expands very rapidly [Hj. 

There is, however, a problem with charmonium suppression as a clear signal of 
a QGP and that is the charmonium suppression observed in pA collisions. In- 
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teractions with co moving particles can break the charmonium apart, the broad- 
ened intrinsic transverse momentum distribution (Cronin effect) and the absorp- 
tion on nucleons also contribute to the suppression. The hard contribution to the 
charmonium-nucleon cross section can be calculated using perturbative QCD and 
is in good agreement with pA data when the formation length and feeding of the 
J / ip from ip' and Xc are taken into account ^3 ^] . A QGP would thus manifest 
itself in an anomalously high charmonium suppression. 

Hadronic cascade models attempt to explain the charmonium suppression in AA 
collisions by interactions with comoving hadrons. They also predict a stronger 
suppression of '0' than oi J /ip [HI • 

The statistical hadronisation model ^H] assumes that all charm quarks are pro- 
duced in hard interactions in the early stage of the collision and are then equili- 
brated in a QGP (thermal but not chemical equilibrium). The charmed hadrons 
are formed at freeze-out according to statistical laws. 

A similar behaviour is expected for the bottonium states. 

Finding signatures that can unambiguously be related to the presence of a QGP 
is very difficult and it is often possible to describe the effects, that where believed 
to be a clear signal, with hadronic scenarios. It is thus likely that the proof for 
QGP formation will be based on several effects, each of which alone cannot provide 
a convincing proof. 

The search at AGS and SPS was based on soft physics such as the strangeness 
enhancement. They are mostly sensitive to the latest stage of the collision after 
hadronisation. At RHIC also hard probes like jet quenching become accessible, 
which provide information on the very early phase of the collision since hard scat- 
terings occur even before the equilibration of the QGP. 
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4 Jet Quenching 



When a hard (i.e. large momentum transfer) scattering takes place in the interaction 
of two protons the two scattered partons leave the protons with high transverse 
momentum. In the centre-of-momentum frame they are emitted back-to-back in 
azimuthal angle due to momentum conservation. The fragmentation of each of 
these energetic partons gives rise to a jet, i.e. a spray of hadrons with small angles 
relative to the momentum of the parton. The energy distribution inside the jet is 
determined by the fragmentation function (Sec. I2.3.1|) . The typical jet signature 
is a large energy deposition localised in a small solid angle. In two-jet events the 
two jets are opposite in azimuth and have nearly the same energy (small deviations 
arise from differences in the momentum fraction x carried by the partons). Partons 
with p± > 3 — 4 GeV are expected to give rise to jets. 

In a soft interaction the cross section for producing a given total transverse 
momentum rises only slowly with the collision energy ^/s since the particle multi- 
plicity rises logarithmically with s and the mean p± also depends only little on s. 
In hard interactions on the other hand the cross section for producing a given total 
transverse momentum rises rapidly with s because the required x ~ ^ \p±\/y/s gets 
smaller and the parton distribution functions rise rapidly towards smaller x. There- 
fore the hard parton-parton scattering dominates the cross section for sufficiently 
high ^ \p±\ and s |20] . 



Figure 4.1: Examples for Feynman diagrams for parton-parton scattering to low- 



The lowest order parton-parton scattering (a + 6 — > c + c?) cross section is of the 
form [20] 





est order 



da 



-{ah — ^ cd) 



1 



\M\' 




dt 
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The quantities with hats refer to the parton level and |A^p is the matrix element 
squared. Spins and colours are averaged in the initial and summed in the final 
state, s, t, u are the Mandelstam variables 



s = {Pa+Pbf = {Pc+Pdf (4.2) 

i={Pa-Pcf = {Pb-Pdf (4.3) 

U= {Pa- Pdf = {Pb - Pcf (4.4) 

Figure ID shows a few examples of processes that contribute to lowest order to the 
parton-parton scattering. 
With the help of the relation 

pl = ^ = LlnH (4.5) 
s 4 

and neglecting the intrinsic transverse momentum of the partons the cross section 
can be written as 

d(T 1 da s da , ^ . 

-jT^ = ^— r = r 4.6 

d{pl) cos^dt t-udt 

Then the lowest order QCD cross section for two-jet production in a collision of 
two hadrons A and B is given by 



%^{AB^ 2 iets) (4.7) 



d(pi) 

1 1 



= j^^"/ '^^^ ® (^^a^fe ~ fa/A{Xa,Q^)fb/B{Xb,Q^)-^^{ab cd) 

dbcd Q Q 

There are many ways to express the differential cross section, a very useful rela- 
tion is 



^d^a d^a d^a , 
E = = 2 (4.8) 

d'^p p±d(j) dy dp± d^dyd {p^ ) 

The agreement between the parton model calculations and data is impressive, a 
few examples are shown in Figure 

Insofar as Bjorken scaling holds in deep inelastic scattering this means that the 
dependence on any dimensionful parameter has disappeared. A similar property in 
hadron-hadron collisions leads to factorised cross sections such that single-particle 
inclusive cross sections are for high pj_ of the form 

E^-^p-^-Fix^,e) x^ = ^ (4.9) 
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where 9 is the centre-of-mass production angle This factorisation in a power 
of and a dimensionless function F is called power law scaling. 

The hadron-hadron cross sections are well described by a parametrisation of the 
form [22] 



E 



dV 



c 



Po 



P±+Po 



(4.10) 



where C, po and n are free parameters. Although it yields a very good description 
of the shape of the cross sections this parametrisation is not suited for extrapolations 
since it doesn't include the dependence on energy and rapidity. 




7t''+X at s'^^=200 GeV 
LO PQCD, K=2/3, <k,">=2.5 GeV' 
NLO PQCD, K=2/3, <k,'>=0.0 GeV' 
NLO PQCD, K=4/3, <k,'>=2,5 GeV' 



2 3 4 5 6 7 8 9 10 11 12 
Pt (GeV/c) 




12 14 
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Figure 4.2: Left: Comparison of pQCD calculations for spectra of charged 
hadrons to data, solid lines are calculations with intrinsic k±, dash- 
dotted lines are without [221 • Right: Invariant cross section for neutral 
pion production, lines are QCD calculations with different parameters 

m 



The hard part of the p_L-spectra in pp collisions is well understood in terms 
of pQCD. Collisions of nuclei are more complicated, nevertheless, because of the 
hard scale involved, high-p^ hadrons are especially suited to investigate nuclear 
modifications. The perturbative processes happen on very short distance and time 
scales so that it is likely that the hard scattering in nucleus-nucleus collisions will 
look exactly like in proton-proton because it doesn't feel the surrounding. However, 
the pdf 's might be different from the ordinary nucleon pdf 's. 

In the following considerations it is assumed that a quark-gluon plasma is formed 
in central nucleus-nucleus collisions. Hard scatterings can thus take place anywhere 
in the overlap region of the colliding nuclei, i.e. inside the plasma. Consequently 
the scattered partons have to traverse it. Since the quark-gluon plasma is a phase 



16 



with very high colour charge density the partons should interact strongly with it. 
The parton, that has a much higher energy than its environment, is expected to 
lose a considerable fraction of it in the plasma. In the subsequent hadronisation 
this leads to a softer jet with lower mean momentum and number of hadrons. 
This expected suppression of high-p^ hadrons as compared to scaled pp yields is 
commonly referred to as jet quenching. 

The amount of energy loss depends on the length of the parton's path trough the 
plasma and is thus given by the geometry of that particular event (i.e. centrality 
of the collision, position of the hard interaction in the plasma, scattering angle of 
the partons). If the energy loss in the plasma is large it should be possible that a 
parton with a long path loses essentially all its energy and does not produce a jet. 
This jet-disappearance would be a clear signal for the QGP, but it is very difficult 
to measure the non-existence of a jet experimentally. A possible way out is offered 
by two-jet correlations. In fact, if the hard scattering took place off the centre 
of the plasma and one parton has a shorter way out so that it loses less energy. 
Consequently the other one has to travel a longer distance through the plasma 
and may get thermalised (i.e. lose essentially all its energy). The parton with the 
shorter path length in the plasma may still have enough energy to produce a jet 
containing high-p^ hadrons. This can then be seen as an event with only one clear 
jet that can be identified experimentally. 

The observation of jet quenching, i.e. suppression of high-p^ hadrons and dis- 
appearance of two-jet events would be a signal for the formation of a QGP. This 
search is currently going on at RHIC (Ch. 16. 2|) . 

A currently accepted scenario for the energy loss mechanism is induced gluon 
radiation. One possible approach to this problem is to consider the scattering off 
static scattering centres [23]. In the approximation that successive scatterings are 
independent and the incident parton has very high energy the radiation spectrum 
induced by multiple scattering can be derived (Sec. 15. 6^ . 

Nevertheless, this may not be the only explanation. It is the aim of this study 
to investigate if soft gluon exchanges can also contribute significantly to the jet 
quenching. A model for this will be dicussed in detail in Chapter ^1 Furthermore 
the general properties of partonic energy loss will be studied in Chapter ^2 The 
main focus will be on the question if the experimental data clearly exhibit signs 
that enable to discriminate between a bremsstrahlung and a scattering dominated 
parton energy loss in the plasma phase. 

It is in fact probable that both effects - scattering and radiation - contribute to 
the overall jet quenching. The parton scatters off a QGP constituent and is deflected 
(and loses energy) which induces gluon bremsstrahlung (and thereby additional 
energy loss). There might also be quantum mechanical interference between the 
two processes. In the gluon bremsstrahlung models the effect of the scattering is 
neglected. The model studied here concentrates in a complementary fashion on the 
scattering and does not take into account gluon radiation. When both effects have 
been studied separately it will hopefully be possible to combine the two into a more 
complete model. 
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5 Basics of Heavy Ion Physics 



Collisions of relativistic heavy nuclei are very complex processes with many inter- 
esting effects that are not well understood today. This chapter cannot cover all of 
these exciting phenomena, instead only the effects that are relevant for this study 
will be presented briefly. 



5.1 Glauber models 

In Glauber models it is assumed that the interactions of two nuclei can be viewed as 
an incoherent superposition of nucleon-nucleon interactions. A collision of two nu- 
clei is then characterised by the mean number of binary nucleon-nucleon collisions 
(iVcou) and the mean number of participating nucleons (A^part)- Both numbers de- 
pend on the impact parameter b. Participating nucleons are those that encountered 
at least one binary collision. 

There exist many variations on the theme, the model introduced in 12^1 will be 
dicussed briefly clS 8b representative. 

The nuclei A and B are characterised by the spherical nuclear density nA.B{i^) 
satisfying 

d^rnA{r)=A (5.1) 
The nuclear density is described rather well by a Woods-Saxon distribution 

^Tir) = , , with no = — -37-^4^— (5.2) 



A simpler alternative that unfortunately doesn't yield as good results is the sharp 
sphere 

<W = 4^0(^^-r^) (5.3) 

The ^-direction is chosen to be the beam axis and the impact parameter b lies 
in the plane perpendicular to z and points from beam to target. The thickness 
function TA{b) is given by 



00 

T4(6) = J dzUAiVWr^) 



(5.4) 
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with 



d%TA{b) = A (5.5) 

Ta can be interpreted as the part of the nucleus A "seen" by a nucleon passing 
through it with impact parameter b. For central collisions (i.e. small values of b) 
Ta is proportional to A^^'^. When going to nucleus- nucleus collisions the product 
of the two thickness functions has to be integrated in the overlap region in order to 
get the nuclear overlap TAsib): 

TAB{b) = jd%d%6\b-h-b2)TA{b,)TB{b2) (5.6) 

with 

Jd%TAB{b) = AB (5.7) 

The product a^^TAsib), where aj^^ is the total inelastic nucleon- nucleon cross 
section at the respective collision energy, can be interpreted as the mean number 
of binary collisions (Neon) at impact parameter b. The total nucleus-nucleus cross 
section is obtained by integrating a^iTAB^b) over b: 

c^i^n^ = / d'baZTABib) = ABaZ (5-8) 

Similarly TAsib) can be multiplied with any cross section to get the mean number 
of events of a particular kind per AB collision at impact parameter b. 
The mean number of participants is given by 

(iVpart(&)) = j d\d%S(b- h - h) {r^(6l)4?(&2) +TB(62)pi1^(6l)} (5.9) 

where p^~^^{bi) is the probability for a nucleon passing through the nucleus A with 
impact parameter bi to take part in at least one interaction. It can be obtained 
from the Binomial distribution: 



= I fA\ f aZTA{b,) y aZTAjb^ 



\ A \ A 



A-O 



1 - ^SS^Mj (5.10) 

~ 1 _ e~<^i^^(''^) for large A (5.11) 

For the most central and symmetric (i.e. A = B) collisions (A^part) is linear in A 
whereas (A'^cou) is proportional to A^^^. It turns out that soft particle production 
in heavy ion collisions scales with the number of participants while hard processes 
scale with the number of binary collisions. 



19 



5.2 Quantum Statistics of an Ideal Quark-Gluon Gas 



As a consequence of the energy dependence of the strong couphng constant 
the interactions among quarks and gluons become weak at high energy densities. 
Thus a gas of non-interacting ultra-relativistic quarks and gluons can be used as 
an approximation. 

The particle numbers are not fixed so that the grand canonical formalism has to 
be used (see e.g. jUj). Let i7 be a set of quantum numbers characterising a one- 
particle state, n(z7) the occupation number of that state and E{r') the corresponding 
one-particle energy. {n{i7)} is the complete set of occupation numbers. Then the 
grand canonical partition function for bosons is 



S(^)(T,V^,/.) = Yl exp(-/3(Etot - /iiV)) 

all states 

{n{u)} \ V 

{n(u)} u 
oo 

= n E - ^,))r^) 

where (3 = 1/T, fi is the chemical potential and Etot and denote the total energy 
and number of particles in the ensemble. The partition function for fermions can 
be derived in a similar fashion. 

1 

u n(i')=0 



Ill + exv(-d(E(u)-u)) ^^'^^^ 



+ exp(-/5(E(i7)-^)) 

lnS(Ty,/i) = ±^ln (1 ± e-^(^(^^)-'^)) (5.14) 

where the plus sign apphes to fermions and the minus sign to bosons. If the energy 

levels are not discrete but continuous the sum can be transformed into an integral 
yielding 

\nE{Ty,^) = ±j In (l ± e'^^^^^^^-^^) (5.15) 
where p and q are momentum and spacial coordinates. 
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The mean occupation number of a state is given by 



d\nE 1 

For the quark gluon gas it is helpful to assume that the QGP is net baryon-free 
as will be the case in ultra-relativistic ion-ion collisions. Then the gluon as well as 
the quark chemical potential vanishes and the numbers of quarks and antiquarks 
are equal. Furthermore, it has to be taken into account that the energy levels are 
degenerate so that the mean occupation number is 

ftp = ^E/T ^ I ^^^'^ ^ ^ \/p'^ + m'^ (5.17) 

where again the plus (minus) sign applies to fermions (bosons) and g is the degen- 
eracy. It is given by the number of flavours times the number of colours times the 
number of polarisations: 



Qg = 2 (polarisation) ■ 8(colour) = 16 (5.18) 
9q = gq = A^/(flavour) ■ 2 (polarisation) ■ 3(colour) = QNf (5.19) 

If the masses are neglected the particle and energy densities can easily be calcu- 
lated: 



^1 = 77 NO / d gd pn„ = / dp- 



V (27r)3 y - ^ - ^ (27r)3 J _ i 



oo 



3a I = _^t32 C(3) ~ 1.2 (5.20) 



27r2 7 ^ e« - 1 27r2 ^ ' ' ' yr^ 



oo 

rig = rig = —j—r:T d qd prip = —— / dp 



V {27^)^ J ^ "^'"^ (27r)3 7 eP/^ + 1 



^T=^c/(3) ~0.9^T=^ (5.21) 



with 



= F(2^F /"'^d'ppn, = ^T'C(4) = (5,22) 
= 5^ T'd(4) = :^r' (5,23) 



n^) = / da — and diO = 77777 / da — (5.24) 

r(0 y e'^-l r(0 y e'^ + l ^ ^ 
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Thus the energy per particle is given by 



^ = 3TM.2.7r (5.25) 

ii = 3T||.3.2T (5.26) 
and the total pressure can be obtained from the equation of state: 

= 1''^^ = + '^ + = ^ + Y^f) T' (5.27) 

5.3 The Bjorken Model 

The Bjorken model is applicable to energetic ion-ion collisions with ^/s > 
100 A GeV. It is based on the assumption that the inclusive particle multiplicity as 
a function of rapidity shows a "central plateau", i.e. particle production looks the 
same in all centre-of-mass-like frames. The second assumption is that the nuclei are 
essentially transparent to each other. The net baryon number is thus contained in 
the receding remnants of the nuclei and the central rapidity region is net baryon-free 
(i.e. the baryon chemical potential vanishes). 

After the collision the two Lorentz contracted nuclei recede from the collision 
point with nearly speed of light. Between them is a dense and hot system that 
expands longitudinally. The expansion velocity is z/t where z = and t = refer 
to the collision point and time. From the charged particle multiplicity observed in 
nucleon-nucleus collisions the initial energy density is estimated: 

. , GeV 
eo = e(ro) ~ 1 " 10 (5-28) 

for an initial proper time 



To = 1 fm/c T = Vt^ - (5.29) 

This relatively high number suggests that the system rapidly comes into local ther- 
mal equilibrium and afterwards follows a hydrodynamic evolution. It is thus possi- 
ble to define a local energy density e(r,y), pressure p{T,y), temperature Tij.y) and 
four- velocity u^ij^y). 

The first assumption implies that the initial condition possesses a symmetry 
under Lorentz transformations which is preserved during the subsequent evolution. 
This means that the energy density, pressure etc. cannot depend on the rapidity. 

Neglecting the viscosity and heat conduction the energy- momentum tensor be- 
comes 

T^^u = (e + p)Uf,u^ - gf,uP (5.30) 
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nonvanishing 
baryon number 



hydrodynamic phase 



Figure 5.1: Space-time diagram of longitudinal evolution of the quark-gluon 
plasma [12] 



It is conserved, i.e. 



dT, 







(5.31) 



Making use of the fact that the occuring quantities depend only on the proper 
time this expression simplifies to 



de 
d7 



e + p 



This can also be written as 

de de dp dT e + p 
dr dp dT dr r 

where s is the entropy density. With the help of 



Ts 

T 



(5.32) 



(5.33) 



dp S ^ de 1 

dT V dp v1 



(5.34) 



{vs is the velocity of sound) the differential equation for the time dependence of the 
temperature is obtained: 



^dT 
T d7 



s 
T 



(5.35) 



In the case of an ideal ultra-relativistic gas the equation of state is e = 3j>. 
Inserting this into Equations 15.321 15.341 and 15.351 leads to the result that the time 
dependence of the energy density and temperature is given by 



e r oc r 



-4/3 



and 

respectively. 



Tir) oc 



(5.36) 

(5.37) 
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The simple picture of purely longitudinal expansion has to be modified when 
the distance between the nuclei becomes comparable to their diameter. There is 
a rarefaction front moving inward from the periphery at the velocity of sound. 
Under the assumption of a time independent velocity of sound the equation for the 
rarefaction front is 



pit) = R 



dt' 



R-VsVt^ - 



(5.3^ 



At transverse distances smaller than the rarefaction front the system continues to 
expand longitudinally because the information that there is a boundary has not yet 
reached this region. Outside the front the gas expands radially outward and cooles 
faster than during longitudinal expansion. 



i ~- ' prompt pions 



nuclear 
"pancakes' 




outward moving, 

more dilute, cooler fluid 



rarefaction front 



Figure 5.2: Geometry of fluid expansion near the edge of the nuclei JH] 



5.4 Azimuthal Anisotropy 

The beam axis [z axis) and the impact parameter h deflne the reaction plane. In 
non-central collisions the overlap region has an almond like shape (Fig. I5.3|l so 
that the pressure gradient depends on the direction. This leads to an asymmetric 
particle emission that can be described by 



E 



d^iV 



d'^p ^■K'pxd'pxdy 



n=l 



(5.39) 



where is the azimuthal angle with respect to the reaction plane [Hj. 

An important feature is that the odd Fourier coefficients change sign at mid- 
rapidity whereas the even are symmetric. It is common to discuss collective flow 
in terms of the lowest order coefficients in the Fourier expansion in Equation 15.391 
They can be determined experimentally in the following way: vi is called directed 
flow and is associated with the mean transverse momentum in the reaction plane: 
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Figure 5.3: Overlap region in a non-central collision 



- ^ (5.40) 
V2 is called elliptic flow and can be written as 

U2 = (cos20) (5.41) 

For negative V2 the emission is dominantly perpendicular to the reaction plane 
whereas a positive V2 favours emissions in the reaction plane. V2 depends on the 
beam energy: It is negative at low energies (BEVALAC/SIS) and becomes positive 
at AGS and SPS with a larger value at SPS. This behaviour is predicted by theory. 
The pressure gradient leads to an emission that is predominantly in the reaction 
plane, but shadowing by spectator nucleons is expected to turn it into an out-of- 
plane emission. This effect becomes unimportant at higher beam energies because 
the spectators are then too far away from mid-rapidity to affect the distributions 
there 



5.5 Cronin effect 

Already in 1975 it was discovered that the production of hadrons with high trans- 
verse momentum is enhanced in proton-nucleus collisions. This effect is called the 
Cronin effect [20] after its discoverer and is commonly quantifled by the Cronin 
ratio R of the p_L-spectra obtained in proton-nucleus collisions using two different 
nuclei A and B {A> B): 

It has roughly the shape shown in Figure 15.41 with a maximum at medium p± . The 
details depend on A, B and the collision energy. 
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Pi Pmax 



Figure 5.4: Approximate beliaviour of tlie Cronin ratio R (Eg. 15.4^ . pi is usually 
of the order 1 GeV and Pmax a few GeV 



The enhancement can be accounted for in terms of multiple scattering of the 
proton as it traverses the nucleus. There are different models jH^ to describe the 
rescattering, but they all lead in one way or another to a broadening of the p±- 
distribution. They differ in the object that experiences rescattering (the proton 
or its partons) and the hardness, i.e. the momentum transfer, of the scattering 
processes. 

A short review of the model introduced in [22] and [SBl will be given here. In 
this case it is the nucleons that undergo rescattering with a predominantly low 
momentum transfer. 

Already in pp collisions it has to be taken into account that the partons have 
an intrinsic transverse momentum kj_ that reflects the size of the proton via the 
uncertainty principle and the Fermi motion. In initial state parton showers further 
transverse momentum can be built up. These two effects are described together in 
a phenomenological approach in EB] • The parton distributions are assumed to 
factorise in two parts that depend on the longitudinal and transverse momentum 
respectively: 



gN{k±,Q'^)fa/N{x,Q'^)dx d^/cj 



(5.43) 



fa/N{x,Q'^) are the normal parton distribution functions [x is the longitudinal 
momentum fraction carried by the parton and Q is the momentum transfer in the 
hard scattering). The transverse momentum distribution gN{k±,Q'^) is assumed to 
have a Gaussian form: 



gN{k±,Q'^) 



1 



7l{kl)N 



(5.44) 



The variance {k^N depends on Q, because both the intrinsic k± and the effect 
of parton showers are included here. It is taken to be 

{kl)N{Q^) = 1.2 GeV^ + 0.2a,(Q2)g2 ^545) 
and the parameters are chosen such that the experimental data are reproduced. 
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In proton-nucleus collisions the proton may experience several soft scatterings 
before the hard process. It is assumed in this model that the k±- distribution is 
still given by a Gaussian. 



In [22] ^a{b) is taken to be {Ncon){b) while it is argued in that faib) cannot 
become larger than 4 because of the dissociation of the proton. But a highly excited 
proton that does not interact as a whole any more will not cease to interact. The 
condition z/^ < 4 is therefore not sensible. The values of 6^ needed to describe the 
data are considerably larger with the restricted Ua, which is a strong indication that 
also the interactions of excited or even dissociated protons have to be taken into 
account. 

The best fit to experimental data is obtained in [22] with 



In nucleus-nucleus collisions the fcx-distributions of both partons are broadened 
accoring to Equation 15.471 so this part is a straightforward generalisation of the 
above results. However, there are more effects that influence the production of 
high-px particles and lead to large uncertainties. An example is the modiflcation 
of the parton distributions in a nuclear environment. 

The Cronin effect is expected to be relatively small in pAu collisions at RHIC 
energies. The different models predict -Rmax — 1.1 ... 1.6 and Pmax — 2.5 ... 3.5 
|31j . Predictions for AuAu are much more difficult due to other nuclear effects as 
mentioned above. 

5.6 Radiative Energy Loss 

A high energy parton traversing a colour charged medium is expected to radiate 
gluons. This medium induced gluon radiation is the QCD analogue to the QED 
bremsstrahlung. There is no unique way to treat this highly involved problem, here 
only a review of the main features of the model introduced in [221 will be given 
without detailed calculations. 

The matter is represented by static scattering centres each of which creates a 
screened Coulomb potential 




(5.46) 



But now the variance is larger due to the soft scatterings: 



{ki)Am = {ki)Nm + - 1) 



(5.47) 






(5.49) 
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where g is the QCD couphng constant and /i is the Debye mass, i.e. the inverse of the 
range of the screened potentiaL It is assumed that the range of the potential is small 
compared with the mean free path A of the parton, i.e. the successive scatterings 
are independent. The calculation of the radiation amplitudes is performed in time- 
ordered pertubation theory. The scattering centres are very heavy so that the 
energy loss in the scattering vanishes. 

Furthermore all partons are assumed to be quarks of very high energy so that 
the soft gluon approximation 




Figure 5.5: Radiation amplitude induced by a single scattering; A, A', B, B', a 
and b denote the colour indices [22] 



The three diagrams that contribute to the gluon emission amplitude induced by 
one scattering are depicted in Figure 15.51 The gluon energy spectrum, which is 
given by the ratio between the radiation and the elastic cross section, is found to 
be 

u-^ = N^^{J{k,qr) (5.51) 
auj 

where J is the emission current 



k± k± - q± 



Jik^<l) = if - : (5-52) 

and 



{J{k,qf) ^ [ d^gl / ,,A k,q? (5.53) 

This result can be generalised to an arbitrary number of scatterings. Figure 15.61 
shows as an example the case of two scatterings. 

The end result becomes quite complicated for a finite number of scatterings and 
the heuristic discussion is more instructive. Here it is assumed that the number of 
scatterings is large 
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Figure 5.6: Gluon emission amplitude induced by two scatterings; a, b and c de- 
note the colour indices 



1 Au2 ^ ^ 

— (5.54) 

There are three different regimes: the incoherent Bethe-Heitler (BH) regime of 
small gluon energies, the coherent regime for intermediate u and the factorisation 
limit corresponding to the highest energies. In the BH regime the radiation is due 
to N incoherent scatterings, in the coherent regime 1 / ^/K scatterings act together 
as a single one (LPM^ effect) and in the factorisation limit finally all centres 
behave as a single scattering centre. This has to do with the formation length of 
the radiated gluons which increases with the gluon energy. The longer the formation 
length is the more scatterings occur during the radiation of the respective gluon. 
The condition 15.541 can also be expressed in terms of the gluon energy: 

ujbu ^ Xfi^ uj u;fact ~ — ^ < E (5.55) 

The last part cufact < E is satisfied when the length L of the traversed medium is 
smaller than the critical length 




L<L,, = J— (5.56) 



^Landau Pomcranchuk Migdal 



29 



Then the radiation spectrum per unit length is 



, UJ < lubu 

" ^BH < ^ < ..fact (5.57) 



L 

The spectrum can be integrated in the range < u < E and < ;z < L in order 
to get the total energy loss. It is found that the energy loss 

AE(L)~a.^(l + o(l)) (5.58) 

is proportional to and independent of E in the high energy limit. In addition 
there is a factorisation contribution which is proportional to asE. For the case that 
L > Lcr or equivalently E < E^r = L'^fJ''^/ A the energy loss becomes 



^E ^as\j^L (5.59) 

Now AE{L) grows only linearly with L but depends on E. The different scenarios 
are reflected in different E- and L-dependencies of the overall energy loss. 
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6 Experimental Results 



6.1 AGS and SPS 

Strangeness Enhancement 

A clear enhancement of strangeness is observed in various systems at AGS and 
SPS. This does not automatically mean that a QGP is formed since the strangeness 
enhancement can partly result from hadronic interactions [S|. 

At AGS an increase of the kaon to pion ratio K/tt is seen and /tx^ > A'^/tt^. 
The data are described by purely hadronic models and the strange particle ratios 
can be fitted by a statistical model with an equilibrated hadronic fireball. This 
does not mean that the system always was in a hadronic phase [3 |H1 Elj . 

At SPS a global enhancement factor 2 for kaons in central PbPb collisions 
is observed and again K~^/n~^ > K~/7r~. The enhancement increases with the 
strangeness content of the respective hadron species [ZlEl- The antihyperon ratio 
H/A is found to smoothly increase from pp over pA to AA collisions. Furthermore a 
strong enhancement of multistrange hadrons from light to heavy nuclei is observed, 
which can be counted as evidence for a non-hadronic enhancement [7j- Hadronic 
transport models fail to describe the SPS data unless they invoke non-hadronic sce- 
narios while microscopic transport models are in good agreement with data [3 |H] . 
The strange particle ratios can be fitted with a hadronic hadron gas and a QGP 

HI IE]. 

This may be seen as evidence for QGP formation, but it has not been possible 
to unambiguously relate the data on strangeness enhancement to the presence of 
a QGP. A more detailed discussion of practical and conceptual problems can be 
found in [7]. 

Chemical Equilibrium and Freeze-Out 

At top AGS and higher energies mid-rapidity particle production can be well de- 
scribed in the framework of a statistical model in full equilibrium while in the 
integrated results strangeness is undersaturated. The freeze-out points are below 
the calculated phase boundary for all except the top SPS energy [1(1] . 

Charmonium Suppression 

J/ip suppression can be quantified using the ratio of the J/ifj cross section to the 
cross section for the Drell-Yan process. The Drell-Yan process is the creation of 
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a lepton pair from qq annihilation {qq It is a hard process that scales 

with the number of binary nucleon-nucleon collisions in a nucleus-nucleus collision. 
The ratio o-jj^/aDY can be interpreted as the probability to produce a. J/tp per 
binary collision. A big advantage is that many experimental uncertainties drop out 
in the ratio. It is often given as a function of L, the mean length of nuclear matter 
traversed by a cc pair ^Hj- It is, however, advisable to be careful since this is not 
a measured quantity and model dependent [7j. 

<^j/i>/ <^dy{L) has been measured at CERN for different systems. The results 
for pA, SU and peripheral PbPb (L < 8fm) are consistent with the exponential 
absorption in nuclear environment while in central PbPb (L > 8 fm) a significantly 
larger suppression is seen [121 El Ell- Hadronic transport models reproduce the data 
but depend heavily on details in the models and need a very high comover density 
[71 12H] which is not typical for hadronic matter. 

The cross section ratio (Tji^/(Jdy{E±) for PbPb collisions as a function of trans- 
verse energy shows a clear suppression in addition to the extrapolated absoption 
in nuclear matter. For peripheral events the data are consistent with the suppres- 
sion expected from normal absorption, but for mid-central collisions the measured 
cross section ratio drops significantly below the normal suppression. It continues 
to decrease for more central events |34j . 

Additional observations are: 

• The loosely bound ip' is similarly absorbed in SU and PbPb [71 115j. 

• The statistical hadronisation model describes the centrality dependence of 
J/ip at SPS, but only with an increased charm production cross section [TH] . 

• The anomalous J/ip suppression observed at SPS can be counted cLS cL si gnal 
for QGP formation, but there are still many uncertainties and the theoretical 
debate is far from settled [3 UHl . 

6.2 RHIC 

At RHIC gold nuclei collide with energies up to 200 GeV per nucleon pair (which 
corresponds to a total c.m.s. energy of ^/s ~ 40TeV). The energy density reached 
is expected to be at least a few GeV / fm^ and thus high enough for the formation of 
a quark-gluon plasma. These calculations are, however, heavily model dependent. 

At these centre-of-mass energies, hadrons with high transverse momentum are 
produced in the fragmentation of partons that underwent a hard scattering. The 
study of these high-p^ hadrons is one of the main goals of the four experiments 
PHENIX, STAR, BRAHMS and PHOBOS. 

The particle multiplicities are so high in ion-ion collisions that jets cannot be 
identified with the normal algorithms which search for a large energy deposition in 
neighboring calorimeter cells. Instead only single high-p^ hadrons can be identified 
and from the relative angle it can be seen whether they belong to the same jet. But 
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since there is a huge background of low-p^ particles it is impossible to identify also 
the softer part of the jet and therefore also the total jet energy as well as the jet 
cone is unknown. 

There have also been p+p and d+Au runs at the same energy (^snn = 200 GeV) 
and a Au+Au run at a/snn = 130 GeV to which the results can be compared. 



6.2.1 Suppression of high-j9^ Hadrons 

In pp collisions hadrons with p± ^ 2 GeV are produced in hard scattering pro- 
cesses. If this is also the case in AA collisions the yield of hadrons with large p± 
is expected to scale with the number of binary collisions (Ch. 15. Ij) . It is therefore 
often characterised by the ratio -Rab of the yield per nucleon-nucleon collision in 
AB collisions to the yield in pp collisions: 

1 d^N^^\ /(A^cou) dVPP 



The inelastic pp cross section is from measurements known to be af^^i ^ 42 mb 
at ^/s = 200 GeV. -Rab is called the nuclear modification factor since all nuclear 
effects will lead to a deviation from unity. (A'cou) has to be calculated for each 
centrality class (i.e. range in impact parameter) separately. It should be noted that 
in the region below ~ 2 GeV the p_|_-spectrum is dominated by hadrons that 
were produced in soft interactions. This contribution is naively expected to scale 
with the number of participants rather than the number of binary collisions. This 
means that even without any nuclear effects -Rab will drop below unity for small 
p±. In summary, the nuclear modification factor is a measure of all kinds of nuclear 
effects and deviations from scaling with (Ncoii) provided hard scattering processes 
are the by far dominant source of high-p^ hadrons. 

The nuclear modification factor is measured by all four experiments jSHl EHl EH 
I2H] • There is no need to discuss all of them since the results are similar and therefore 
only the PHENIX data will be presented here. PHENIX measures charged hadrons 
and neutral pions covering the pseudorapidity range |?7| < 0.35, but for the charged 
particles a cut of |?7| < 0.18 was applied. The events are divided into nine centrality 
classes that are listed in Table IHUl together with the corresponding values of (A^'cou) 
and (A^part)- 

The charged particle spectrum measured by PHENIX in the pp run does not 
reach far enough in p±, so the pp reference spectrum was constructed from the vr" 
spectrum measured by PHENIX and the observed charged hadron to neutral pion 
ratio (PHENIX and other experiments ^^)- 

The results for -RauAu at ^snn = 200 GeV are shown in Figure 16.11 There is a 
clear suppression of high-p^ hadrons that gradually increases with increasing cen- 
trality. The event centrality is determined from the number of spectator neutrons 
and the number of fast secondaries. The drop in the region p± ^2 GeV is certainly 
caused at least partly by the (Acoii) instead of (A^part) scaling in Equation 16. II The 
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Table 6.1: Centrality classes with the number of participants and binary collisions 
as used by PHENIX, the percent values characterising the centrality 
classes give the fraction of the total cross section ,7^"^" = 6.9 b 23 



centrality 


(iVcou) 


(iVpart) 


0- 


10% 


955.4 ±93.6 


325.2 ±3.3 


10- 


20% 


602.6 ±59.3 


234.6 ±4.7 


20- 


30% 


373.8 ±39.6 


166.6 ±5.4 


30- 


40% 


219.8 ±22.6 


114.2 ±4.4 


40- 


50% 


120.3 ± 13.7 


74.4 ± 3.8 


50- 


60% 


61.0 ±9.9 


45.5 ±3.3 


60- 


70% 


28.5 ±7.6 


25.7 ±3.8 


70- 


80% 


12.4 ±4.2 


13.4 ±3.0 


80- 


92% 


4.9 ± 1.2 


6.3 ± 1.2 


min. 


bias 


257.8 ±25.4 


109.1 ±4.1 



ratio is flat for pj_ > 4.5 GeV indicating that particle production is dominated by 
hard scattering in this regime although the yield per binary collision is not as high 
as in pp collisions. The scaling with x± {x± = 2p±/^/s) gives further support to 
this hypothesis jHS]- The increase of the suppression with centrality is consistent 
with a plasma scenario since the size and thereby the energy loss suffered by hard 
partons are also expected to increase with centrality. From the Cronin effect an 
enhancement of particles with intermediate p± (Ch. 15. 5p is expected so that the 
suppression is possibly even higher. 

It has been argued in [221 that in case of plasma formation and high energy loss 
in the deconfined medium only the partons produced near the surface could escape 
since they would be thermalised if the way through the plasma was too long. This 
would result in a scaling with (A^'part) instead of {Ncow) also for the high-p_L hadrons 
|39j . There is, however, no clear evidence for that. 

A surprising result is the centrality dependence of the nuclear modification factor: 
It scales neither with (Ncoii) nor with (A^part) but is linear as a function of percent 
of the total cross section. Figure 16.21 shows the centrality dependence of the mean 
-RauAu for P± > 3 GeV of the PHENIX vr*^ data, the linear fit is 

(i?AuAu(p± > 3GeV)) = (0.00835 ±0.00034) -centrality [%] ± (0.186±0.012) (6.2) 

There is no convincing explanation for this behaviour so far. 

An astonishing feature of Figure 16.11 is the stronger suppression of neutral pions 
at intermediate p±. It turns out that this is not a stronger suppression of vr" but 
a higher production rate of protons and antiprotons. Figure 16.31 shows the ratio of 
charged hadrons to neutral pions in AuAu for central and peripheral collisions. The 
dotted hne at 1.6 indicates the ratio measured in e~^e~ and pp collisions. This plot 
nicely illustrates the excess of charged particles at intermediate p±. From Figure 
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Figure 6.1: RauAu{p±) for vr'' and charged hadrons at mid-rapidity for the different 
centrahty classes measured by PHENIX 



it becomes clear that this is caused by an enhanced (anti)proton production in more 
central AuAu collisions indicating a deviation from the standard picture of particle 
production in the medium-p^ range (2 GeV < p± ^ 4.5 GeV). 

In deuteron-gold collisions at RHIC no suppression of high-p^ hadrons is ob- 
served. At mid-rapidity there is an enhancement [HTl EUJ |^ that is usually in- 
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Figure 6.2: Centrality dependence of the mean -RauAu for p± > 3 GeV for 
PHENIX 7r° data [SI] 
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Figure 6.3: Charged hadron to neutral pion ratio in central (0 — 10% - squares) 
and peripheral (80 — 92% - circles) AuAu collisions (the peripheral 
data points are offset by 130 Me V); PHENIX 



terpreted as the Cronin effect while the ratio is consistent with unity at larger 
rapidities jlSj- This further strengthens the case for a plasma scenario since the 
suppression in AuAu collisions is obviously related to final state effects (initial state 
effects should be present in the dAu data as well). Figure shows the nuclear 
modification factor -RdAu in dAu collisions measured by PHENIX, the values of 
(-/Vcou) and (A^part) are 8.5 and 1.7 respectively. 

In summary, there is a substantial suppression of high-p^ hadrons in AuAu col- 
lisions that smoothly increases with increasing centrality. The spectral shape is 
unchanged as compared to pp for large p± indicating that fragmentation of hard 
partons is the dominant particle production mechanism in this regime. At medium 
p± a large relative proton and antiproton yield is observed, which is a hint to a 
deviation from the "normal" particle production. All observables approach their 
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Figure 6.4: p/ir (left) and p/ir (right) ratios; open (filled) points are for charged 
(neutral) pions; PHENIX ^ 
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Figure 6.5: -RdAu for charged hadrons and neutral pions at a/snn = 200 GeV; 
PHENIX HOj 



pp values in the most peripheral events. It should be kept in mind that the scaling 
behaviour is different for small p±. No suppression is observed in dAu collisions 
indicating final state effects as the most probable origin. These results are consis- 
tent with a scenario where a quark-gluon plasma is formed in the overlap region of 
the colliding nuclei and energetic partons lose a substantial amount of energy when 
traversing the plasma. However, this is an indication and no proof. 

Surprisingly the nuclear modification factor increases linearly with centrality, this 
feature is presently not understood. 
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6.2.2 Disappearance of Back-to-Back Hadron Correlations 



STAR also measures the azimuthal correlation of charged particles at mid-rapidity 
(|?7| < 0.7) jmiH]. The azimuthal distribution is defined as 

^^^^^ = /V^H^ ^6-3) 

Particles with 4 < < 6 GeV are defined as trigger particles. For each trigger 
particle the azimuthal separation A0 from all other particles satisfying 2 GeV < 
P± < P±(trig) is calculated yielding iV(A0). In pp collisions D{A(f)) shows clear 
jet-like peaks at A0 ~ and A0 ~ tt (Fig. 16. (i1 and IFTTT a)). The distribution is 
characterised by a Gaussian for each peak: 

where the indices and B stand for the near-side and back-to-back respectively 
and P is a constant offset. The fit is shown as the solid line in Figure I^TT a). the 
parameters are given in Table | 



Table 6.2: Fit parameters from Equation 16.41 for the AuAu central sample only 
one Gaussian was fitted to the data after subtraction of the fiow and 
pedestal contribution [l^ 





pp min. bias 


dAu min. bias 


AuAu central 


An 


0.081 ±0.005 


0.074 ±0.003 


0.093 ±0.008 




0.18 ±0.01 


0.20 ±0.01 


0.22 ±0.02 


Ab 


0.119 ±0.007 


0.097 ±0.004 


/ 




0.45 ±0.03 


0.48 ±0.02 


/ 


p 


0.008 ±0.001 


0.039 ±0.001 


0.004 ±0.003 



In nucleus-nucleus collisions there will also be a contribution from the elliptic 
fiow anistropy (Ch. 15. 4j) : 

|^ = fi[l±2.;^os(2A0)] (6.5) 

where f 2 has been measured independently in the same set of events and is assumed 
to be constant in 2 < j)^ < 6 GeV, whereas B has to be fitted. The values for the 
different centrality classes are listed in Table 

Figure 16.61 shows the measured azimuthal distribution for different centralities 
in AuAu together with the elhptic flow contribution (solid lines) and the pp data. 
Jet-like correlations can be seen as an excess of the measured correlation over the 
elliptic flow contribution. The jet on the trigger side (A0 ~ 0) is present for 
all centralities but the away-side jet (A0 ~ vr) vanishes in more central events. In 
Figure|niZ(b) the distribution in pp collisions is compared to central AuAu collisions 
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Table 6.3: Values of A'part, V2 {2 < p± < 6 GeV) and the normalisation constant 
B for the different centrality classes 



centrality 


-^part 


V2 


B 


0-5% 


352 ±7 


0.07 ±0.01 


1.442 ± 0.003 


5-10% 


298 ± 10 


0.10 ±0.01 


1.187 ±0.008 


10-20% 


232 ± 11 


0.15 ±0.01 


0.931 ±0.006 


20-30% 


165 ± 13 


0.19 ±0.01 


0.633 ±0.005 


30-40% 


144 ± 13 


0.21 ±0.01 


0.420 ± 0.005 


40-60% 


61 ± 10 


0.22 ±0.01 


0.231 ±0.003 


60-80% 


20 ±6 


0.24 ±0.04 


0.065 ±0.003 
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Figure 6.6: Azimuthal distribution D{A(j)) for pp and different centralities in 
AuAu measured by STAR, the lines show the elliptic flow and pedestal 
contribution (Eq. 16. 51 Tab. 16. 3|) 



after subtraction of the elliptic flow and pedestal contributions. It is not surprising 
that the correlation on the trigger side is nearly as strong as in pp collisions since 
the trigger condition introduces a bias. The important point is that the jet on the 
other side disappears completely in the most central collisions. The suppression 
increases again smoothly with centrality. 

The results from the dAu measurements are shown in Figure lUlTT a) and Table 
There is no substantial suppression of the away-side jet which again points towards 
final state effects being responsible for the disappearance of the second jet in central 
AuAu collisions. 

The results are commonly interpreted in terms of a plasma scenario in the fol- 
lowing way. The data samples for the AuAu collisions are heavily biased towards 
events where the hard scattering took place near the surface of the overlap region 
due to the trigger condition. Particles with p± > 4 GeV are rare already in pp 
collisions. If the energetic parton is produced inside the plasma and has to traverse 
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Figure 6.7: Azimuthal distribution D^Acf)) for pp, dAu and central AuAu (elliptic 
flow and pedestal distributions have been subtracted in (b)), solid 
lines in (a) are Gaussian flts (Eg. 16. 4^ Tab. 16.2(1 . STAR jH] 



a substantial part of it, it will lose a sizable fraction of its energy in interactions 
with the plasma. The leading hadron formed in the fragmentation will then have a 
transverse momentum which is too small to pass the trigger condition. If now the 
hard scattering happened near the surface with one parton escaping undisturbed 
the other one will typically have to traverse a substantial part of the plasma for 
geometrical reasons. The mean path length through the plasma increases with 
increasing centrality and one would therefore expect an increasing suppression of 
the away-side jet. In dAu collisions, where no plasma is formed, the azimuthal 
correlation should be similar to pp distribution. 

This argumentation sounds convincing, nevertheless one has to be careful. Firstly 
the p_L-spectrum has a tail to very high energies and these partons may produce 
trigger particles even after heaving lost some energy. Therefore the trigger bias 
does not automatically mean that the hard scattering took place near the surface. 
This point can only be clarifled with the number of trigger particles produced per 
binary collision which is much smaller in the scenario where only scatterings that 
happen to be located near the surface can give rise to trigger particles. Secondly 
the lifetime of the QGP which is estimated to be of the order of a few fm limits 
the path lengths. In the surface emission scenario this would lead to saturation as 
soon as the geometrically available path length becomes shorter than the lifetime. 
And flnally the matrix element is symmetric in azimuth which means that there is 
a sizable probability that both partons escape nearly undisturbed when they are 
emitted from the surface. Consequently the observed azimuthal correlation seems 
to contradict the surface-emission scenario and the observed jets stem from some 
greater depth. 
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The conclusion can be drawn that also the measured azimuthal correlation of 
high-px hadrons is qualitatively consistent with the hypothesis of quark-gluon 
plasma being formed in relativistic nucleus-nucleus collisions. 



6.2.3 Correlation with the Reaction Plane 

The almond-like shape of the overlap region in mid-central collisions leads to a 
dependence of the two-particle correlation (Eq. Ifj.Hj) on the orientation of the pair 
relative to the reaction plane. A parton that is emitted in the reaction plane has 
on average to traverse a shorter distance in the QGP than one that is emitted 
perpendicular to it. Consequently the suppression of the away-side jet is stronger 
in the latter case. 

The effects of different emission angles with respect to the reaction plane have 
alse been studied by the STAR experiment [33] . The procedure is analogous to the 
one described in Section f6.2.2l but here all tracks with |?7| < 1 are considered. A 
pair is defined to be in-plane if the azimuthal angle of the trigger particles satisfies 
|0t"g — VI/2I < 7r/4 or 10*"^ — ^1/2! > 37r/4 where is the reaction plane angle of the 
event. A trigger is out-of-plane if 7r/4 < — \l/2| < 37r/4. 

The elliptic flow contribution is [32] 



^ out 



dA0 



trig 
TTV. 



2^assoc ( -'2 
71 



±2 cos 2A* \ , , 

—r^ ^ ^ COS 2A0 

±4t;fs(cos(2A*))y ^ ' 



(6.6) 



where vf^^°'^ and v^^^ are the elliptic flow of the associated and trigger particles, 
respectively, and (cos(2A\l/)) is the reaction plane resolution. The values for 20%- 
60% centrality are = 0.20, vl"^ = 0.18 and (cos(2A^)) = 0.70. 

The azimuthal distributions are shown in the upper panel of Figure l6^ The solid 
symbols are data points whereas the open symbols are reflections around A0 = 
and A(f) = ir. The elliptic flow contributions from Equation 16.61 were fitted with B 
as the only free parameter. The results, shown as dashed lines, are B = 0.649 ± 
0mA{stat.)±0m5{sys.) for the in-plane and B = 0.638±0.004(stat.) ±0.002(s2/s.) 
for the out-of-plane distribution. 

The lower panel of Figure 16.81 shows the distributions after subtraction of the 
elliptic flow together with the pp measurement. The data clearly confirm a stronger 
quenching of jets perpendicular to the reaction plane where the original parton had 
to travel a longer distance inside the QGP and thus experienced larger energy 
loss. The large uncertainty is mainly due to the uncertainty in the elliptic flow 
contribution. 

The fact that a substantial part of the away-side jets survives in the reaction 
plane also illustrates that the QGP cannot be completely opaque since the smallest 
diameter of the overlap region in a 20% centrality event is already ~ 4 fm. It can 
thus be concluded that a considerable part of the jets survives a journey of a few 
fm through the QGP. 
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Figure 6.8: Upper panel: Two-particle correlations for trigger particles in-plane 
(squares) and out-of-plane (triangles) in mid-central AuAu. The con- 
tribution from elliptic flow is shown by dashed lines. Lower panel: 
Azimuthal distributions after subtracting the elliptic flow and pp 
data. I 



6.2.4 Other Observables 

The yields of multi-strange particles are enhanced as compared to SPS. The data 
indicate strangeness saturation for central collisions [IHj. 

The particle ratios are again beautifully described by a statistical model assum- 
ing full chemical equilibrium. The freeze-out point is found to be consistent with 
the phase boundary calculated with lattice QCD jlTj. Furthermore it has been 
shown that the chemical freeze-out temperature coincides with the critical temper- 
ature since it is not possible to maintain chemical equilibrium in a hadronic phase 
even close to the phase transition |1H]. This provides a strong argument for QGP 
formation. 

The first measurements of J ftp production show a suppression but the experi- 
mental uncertainty (mainly lack of statistics) is still too large to make more detailed 
statements |inj. 



42 



7 Comparison of Identified Particle 
Spectra from SPS and RHIC 



While the hard part of the p_L-spectra in nucleon-nucleon colhsions is theoretically 
under control and shows the power-law behaviour typical for hard parton-parton 
scattering this is not the case for the soft part. Here the main contribution comes 
from non-perturbative processes that are not well understood so that one has to 
rely on paramtrisations of measured cross sections. However, the situation is not as 
bad in ion-ion collisions where the soft part of the p_L-spectra is better understood 
and contains valuable information about the properties of the created medium. 

In collisions of relativistic nuclei the \ow-p± part of the hadron spectra is de- 
termined by the hydrodynamic behaviour of the system, primarily the freeze-out 
temperature and the pressure which drives the transverse expansion. A still open 
question is whether the physics at SPS energies is qualitatively different from the 
physics at RHIC. The hard part of the spectra is hardly accessible at SPS, but 
the soft part can be used to quantify differences between RHIC and SPS. In this 
chapter a detailed comparison of the shapes of the p±- spectra of identified particles 
is carried out. 

A thermal source at rest emitting particles from the origin leads to a spectrum 
of the form 

oce T (7.1) 



'm±dm±dy 

which is completely determined by the temperature. m± = \/px + rnl is the trans- 
verse mass and niQ is the respective particle's rest mass. The transverse expansion 
leads to deviations from this form which are, however, model dependent. Further 
distortions arise from resonance decays. 

The freeze-out temperature can be extracted from the ratios of different particle 
yields and is found to be nearly the same for SPS and RHIC, namely T = (168 ± 
2.4) MeV at SPS HOI and T = (174 ± 7) MeV at RHIC HH. Differences in the 
shape of the spectra should thus stem from different expansion velocities. 

In order to avoid the problem with the model dependence of the transverse ex- 
pansion contribution an exponential aexp {—b{m_i_ — mo)} is fitted locally, i.e. in 
a small p_L-interval, to the p_L-spectrum. b will then depend on p± and cannot be 
interpreted as an inverse temperature, but it can be used to quantify differences in 
the slopes of the spectra. Since this comparison relies solely on the shape of the 
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spectrum the parameter a which describes the yields is irrelevant. To get compara- 
ble results the spectra from central collisions near mid-rapidity were used. A large 
number of different identified particle spectra were analysed in order to improve 
the significance and give an overview over the present status in terms of experi- 
mental data. All analysed spectra are listed in Table mi Concerning protons and 
antiprotons only spectra that are corrected for feed-down from week decays were 
considered. Figure 17.11 shows as an example a measured pion spectrum with the 
fitted exponentials. 




1e -14 I 1 1 1 < <- 

2 4 6 8 10 

rriT-mg / GeV 



Figure 7.1: vr*^ spectrum measured by the PHENIX experiment [21] with locally 
fitted exponentials 

A potentially critical point is the possible dependence of the results on the in- 
terval length. Fortunately it was found be small as is illustrated in Figure 17.21 A 
compromise between a small interval, which is desirable because of better sensitivity 
and resolution, and a small error has to be made. The interval length was there- 
fore chosen individually for each spectrum depending mainly on the experimental 
uncertainties. 

The results for the kaons are shown in Figure 17^ as a representative, the rest can 
be found in Appendix^ The fact that b changes with p± illustrates that transverse 
expansion does indeed play a role. The differences between SPS and RHIC are 
small, the RHIC points are continuously slightly higher. The discrepancies among 
the experiments at the same accelarator are also sizable, in fact they are of the 
same magnitude as the differences between RHIC and SPS. There is also no visible 
difference between the SPS data at 80 A GeV and 158 A GeV beam energy, except 
maybe in the A spectra. The sum of the results clearly suggests that there is no 
dramatic change in (soft) physics when going from SPS to RHIC. There is room 
for a small quantitative difference. 
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Figure 7.2: Fit parameters for tiie vr*^ spectrum from Figure 17.11 obtained with 
different interval lengths 



PHENIX 200AGeV K* i — ^ 
PHENIX 200AGeV K" 
STAR 200AGeV K* 
NA57 158AGeV K° 
NA49 158AGeV K+ o 
NA44 ISBAGeV K* » 
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Figure 7.3: Fit parameters for the kaon spectra 



A comparison of the elliptic flow leads to the same conclusion, here the differences 
between SPS and RHIC are also not very big. Indeed the two observables are 
correlated because the pressure, that drives the transverse expansion, also has an 
influence on the pressure gradient, which is responsible for the elliptic flow. The 
transverse size of the system also turned out to be similar at SPS and RHIC [S2]. 
Together these observation lend even stronger support to the conclusion that at 
least the soft physics is similar at SPS and RHIC. 
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Table 7.1: List of all analysed spectra 
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8 Monte Carlo Techniques 



A very general definition would be that Monte Carlo is everything that makes use 
of random numbers to solve a problem. Monte Carlo methods are mainly used to 
simulate processes, estimate integrals or select random variables according to some 
probability density (or a combination of these). The processes to which Monte 
Carlo methods are applicable do not necessarily have to be of stochastic nature, 
although this is normally the case in high energy physics. 

The values that a random variable takes are unpredictable, so there have to be 
at least two different possible values of that variable. The distribution of a random 
variable is given by a probability density g{u)du, which is the probability that the 
value lies in the range [u,u + du]. The scattering angle inaa + 6— i>a + 6 elastic 
scattering is an example for a random variable. In this case the probability density 
is given by the differential cross section da = f{'d)d'd. This example also illustrates 
that probability densities in physical problems are often not normalised to unity, 
because the integral may also have a physical meaning (in this example the total 
cross section). 

The integral of a function f{x) in an interval [a,b] can be estimated by choosing 
values Xi distributed randomly in [a,b], summing the function values f{xi) in 
these points and dividing by A^. According to the law of large numbers the Monte 
Carlo estimate of an integral converges (in the statistical sense) to its exact value 
as the number of points approaches infinity. 



Convergence in the statistical sense means that for every probability p and pos- 
itive number e there is a /c such that for all N > k the probability of the difference 
between the left- and the right-hand side of Equation 18.11 being smaller than e is 
greater than p. 

The Monte Carlo estimate is a function of the random numbers {xj}. For finite 
A^ the value obtained with a different set of numbers {x'j} will be different from the 
first one. The central limit theorem says that for large A^ the Monte Carlo estimate 
is a Gaussian distribution. Its expectation is the true value of the integral. The 
expectation of a function f{x) is defined as 





a 
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and the variance (average of the squared deviation from the expectation) is 




(8.3) 



The variance of the Monte Carlo estimate is V{f)/N. The speed, with which 
the Monte Carlo estimate converges, does not depend explicitly on the number of 
dimensions (although other calculations in the algorithm might take more time in 
higher dimensions). This means that for small number if dimensions numerical 
quadrature will be much more efficient while the Monte Carlo clearly wins for 
higher number of dimensions. Since the variance also depends on the variance 
of the function /, the performance of the Monte Carlo can also be improved by 
reducing V{f). These methods will not be discussed explicitly here but can be 
found for example in jHl]. The variance reducing techniques described below for 
the selection from a given probability density can be adopted to the estimation of 
integrals. 

Given a random number generator that produces numbers R uniformly dis- 
tributed in ]0,1[ these can be used to obtain numbers distributed according to any 
other density function f{x) defined on an interval [a,b]. This is extremely useful in 
simulations of physical processes. 

The simplest way to do that is the hit-or-miss Monte Carlo. First an upper 
limit /max of the function has to be found, i.e. /max > f{x) Vx G [a,b]. Then the 
procedure is to select an x from a fiat distribution, i.e. x = a + R{h — a) where R is 
a random number in ]0,1[. Then another random number R in the interval ]0,1[ is 
selected and if i? < /(a;) //max the x- value is accepted. Otherwise it is rejected and 
the procedure continues with choosing a new x. 

The advantage of this algorithm is that it works for practically all functions 
without singularities. The efficiency is 



If / varies a lot the efficiency becomes very bad, so it might be useful to divide the 
interval into suitable subintervals. Then a decision which subinterval to choose has 
to be made first (the relative probability is given by the integral in the subinterval). 

Another method is the direct sampling, that can be used if / has a primitive 
function F whose inverse can be calculated. Then x-values can be found using 
the relation 



X 



b 




(8.5) 



a 



a 



=>x = F-\F{a) + R{F{b) - F{a))) 



(8.6) 
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This is a very elegant and efficient method with the main drawback that the number 
of functions to which it can be apphed is very hmited. 

The efficiency can also be improved by generating more points in regions where 
the function values are large (importance sampling). This is mathematically equiv- 
alent to a change of integration variable: 



In practice it works like this: First a suitable function g{x) has to be chosen. 
Suitable means that it satisfies (i) g{x) > f{x)\/x G [a,b], (ii) has an invertible 
primitve function or a ^f-distributed random number generator is available and (iii) 
the ratio f{x)/g{x) should be as constant as possible. Then select a x according to 
g{x). Then a random number R is selected, the corresponding x G [a,b] is accepted 
if -R < f{x)/g{x) and rejected otherwise. The importance sampling will improve 
the efficiency if V{f/g) < V{f), but it suffers from the small number of available 
functions g. 

A crucial point in all Monte Carlo calculations or simulations is the random 
number generator, since reliable random numbers are needed. Truly random num- 
bers are unpredictable and therefore also unreproducible. They can be obtained 
from physical processes such as radioactive decay or thermal noise in an electronic 
circuit. However, these truly random numbers have some disadvantages when it 
comes to their application. One point is that it is difficult to construct a device 
that generates random numbers by observing suitable processes being at the same 
time accurate and very fast. In fact speed is of great importance, since extended 
simulations consume huge amounts of random numbers. 

An alternative to using truly random numbers is to generate pseudo-random 
numbers from a mathematical formula. These are not random in the sense that they 
are reproducible, but are just as good for Monte Carlo calculations as long as they 
are uncorrelated (which is of course the critical point). There are also advantages 
in the sequence of numbers generated being always the same for one particular 
starting point, e.g. debugging of a program using random numbers becomes much 
easier. When generated on a computer the sequence of numbers has a finite length 
and at some point there will be a number that has already occured before. From 
now on the sequence will repeat itself and therefore be useless. So the period is an 
important characteristic of a (pseudo-) random number generator. The generator 
used in this study has a period of over 10^^. 

The ideal algorithm for producing pseudo-random numbers does not exist, so 
testing is an issue. But even with a lot of tests made it is never possible to be sure 
that the numbers are suitable for a certain problem that is to be investigated using 
Monte Carlo methods. It should therefore be kept in mind that the generator one 
is using is not perfect and that problems may occur. The lesson to learn is that 
random number generators have to be handled with care. 




where 
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Monte Carlo event generators are powerful tools for the simulation of compli- 
cated processes like particle collisions, where one has to deal with many degrees 
of freedom. Instead of trying to solve the whole problem in one step (which is 
practically impossible), the process is divided into a sequence of subprocesses each 
of which is of reasonable complexity. For high energy physics, for instance, the 
first part consists of the particles before the interaction, where they are treated 
in terms of their parton density functions. The next step is the hard interaction 
that can be calculated using pQCD matrix elements. Then parton showers a la 
DGLAP are added, after that the event is hadronised and as a last step follows the 
decay of unstable particles. In each part the dynamical variables are chosen from 
their probability distributions. The outcome of this whole machinery is a list of all 
particles in the final states with their masses and momenta, but also the complete 
history of the event is provided. When generating many events the cross sections 
can be estimated at the same time. Ideally the mean and variance of a sample of 
generated events should be the same as for the real physical process. 

For this study the Pythia 6.2 package [0] was used, but with an additional 
phase, representing the interaction of the scattered partons with the plasma, added 
between the generation of parton showers and hadronisation. 
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9 Soft Colour Interactions 



As was already discussed in Chapter 12. II only QCD processes involving a large mo- 
mentum transfer can be treated in pertubation theory, since the coupling strength 
as increases with decreasing energy. This is exactly the reason why it cannot be 
assumed that interactions below the cutoff for perturbative calculations do not oc- 
cur. They may be negligible in many situations but there are observables related 
to the hadronic final state that are sensitive to the soft interactions. 

The Soft Colour Interaction (SCI) model P is an attempt to gain a better un- 
derstanding of non-pQCD (non-perturbative QCD) processes. It describes the in- 
teraction of partons originating from a hard scattering in ep or pp collisions with 
the colour background field inside the hadron. The basic assumption is that the 
change in the dynamics caused by the soft interactions is negligible. The important 
point is that the partons can exchange colour leading to a different string topology, 
which may then result in a different final state after hadronisation of the string. 

The background field is represented by the remnants of the hadron(s), i.e. those 
consituents that did not take part in the hard scattering (Ch.|21). In the case where a 
valence quark was kicked out the remainder is a diquark in a colour antitriplet state. 
If the scattered parton is a gluon the remnant (namely the three valence quarks 
in a colour octet state) is split into a quark and a diquark. If a seaquark was 
scattered its partner has to be taken into account to conserve quantum numbers. 
The remnant then consists of the three valence quarks (again split into a quark and 
a diquark) and the seaquark partner. Together they form a (anti)triplet, depending 
on whether the seaquark or the sea-antiquark is left. 

The interactions consist of an exchange of a colour octet (consisting of colour 
and anticolour) between a pair of partons, i.e. this can be viewed as a soft gluon 
exchange. The model is implemented as add-on routines to Pythia at a stage 
between the hard scattering process including parton showers and the hadroniza- 
tion. Each perturbativly produced parton may interact with each of the remnant 
consituents with a probability Pint which cannot be calculated from theory and is 
the only free parameter of the model. The outcome of this procedure will vary from 
event to event. 

One example for a deep inelastic scattering event is shown in Figure 19.11 In the 
conventional Lund model the remnant is in a colour octet state with the colour 
and anticolour charges matching those of the gluon so that gluon and remnant 
together form a singlet. In the case shown in Figure WJ\ the radiated gluon carries 
the anticolour charge of the original gluon (which belongs to the charge of the 
quark in the remnant) and its colour is the anticolour of the original gluon after 
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Figure 9.1: Photon-gluon fusion in deep inelastic scattering with the strings drawn 
as double dashed hues; left: configuration in the Lund model, right: 
configuration after the exchange of a soft gluon (colour octet) between 
the quark from the 2 — > 2 scattering and the quark in the remnant 

in 

the radiation. When this gluon splits into a qq pair the quark gets the colour and 
the antiquark the anticolour charge resulting in the string configuration shown in 
Figure I^TTl The exchange of a soft gluon between the two quarks in the diagram 
leads to an interchange of the two colour charges resulting in a different string 
topology (right-hand side of Figure WA\ . Now the remnant system is in a singlet 
state. If instead the antiquark had interacted with the diquark, the scattered quark- 
antiquark pair would form a colour singlet system and the remnant, together with 
the radiated gluon, would form the other singlet. 

The Generalized Area Law (GAL) model [02] is formulated in a string basis 
instead of a parton basis, since the soft gluons may not resolve single partons. Here 
the interactions are between strings and constructed such that the area swept out 
by the strings in energy-momentum space is minimised. The basic idea, however, 
is the same as in the SCI model. 

The two models give very similar results, although there are differences when 
attention is paid to the details. In this study only the SCI model was used. The 
parameter Pint was chosen to reproduce the rate of diffractive rapidity gap events 
observed at HERA (a rapidity gap is an interval in rapidity or polar angle where 
no particles are produced) [T]. Its value is Pint = 0.5, but the observables studied 
so far do not depend strongly on the exact value. Once Pint is fixed the models 
describe successfully not only the diffractive data from HERA, but also production 
of W , jets and beauty in diffractive events and double leading proton events at the 
Tevatron p. Furhermore it provides a very general description in the sense that 
it does not introduce substantial changes in non- diffractive events (events without 
rapidity gaps). 

The SCI model has recently received a firmer theoretical basis in terms of parton 
rescattering in QCD [HH] . 
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10 The Soft Colour Interaction Jet 
Quenching Model 



10.1 Description of the Model 

The model is built on the same concept of non-perturbative QCD-interactions be- 
tween (colour-) charged objects as the SCI model (Ch. IH)). But now a scattered 
parton has to traverse a part of the QGP instead of only the proton remnant and 
will thus experience much more interactions. It is therefore necessary to consider 
also the momentum-transfer since even a small energy-loss per interaction may sum 
up to give a sizeable amount in the end. 

The simulation is based on proton-proton collisions in Pythia, where the addi- 
tional interactions with the plasma are added after the perturbative part but before 
hadronisation. This is justified since at ^/s = 200 A GeV the probability to have 
more than one hard scattering per ion-ion collision is very small. All soft processes 
like hadronisation of the plasma, collective flow etc. are not included in the model 
which is only concerned with phenomena that are connected to high-p^ partons. 
The only non-perturbative process in the model is the interaction of these partons 
with the QGP. In the hard scattering a transverse momentum transfer of at least 
3 GeV is required. It is necessary to have a p_L-cut since the matrix element is ap- 
plicable only in the perturbative regime and diverges aX p± — 0. Here a relatively 
high value is used because non-perturbative interactions influence the particle pro- 
duction up to ~ 2 GeV. Furthermore the value of 3 GeV allows to accumulate a 
good statistics in the high-p^ regime up to ~ 10 GeV in reasonable run times. 

10.1.1 Geometry 

The coordinate system is chosen such that the origin is in the centre of the overlap 
region, the z direction corresponds to the beam direction and the x direction is 
defined by the impact parameter fPig. 110. T|) . 

The diameters of the overlap region in x and y direction are 

4 = 2i?Au-fe (10.1) 

dy = ^J^Rl^ - (10.2) 

The geometrical effects are modeled with the help of a simple Glauber model 
(Sec. 15. ll [2ni) with a sharp sphere potential. The centrality (defined as fraction of 
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Figure 10.1: Definition of tlie coordinate system 



tlie total cross section) is proportional to since da ochdh for geometric reasons. 
From a fit to a Glauber model calculation [2E] 

h = [(1.579 ± 0.003) fm^ • centrality [%]] (10.3) 

is obtained. In the model a centrality range has to be given which is then translated 
into a range of the impact parameter and sampled according to dcr oc bdb. 

For the Monte Carlo simualtion it is necessary to have a simple parametrisation 
of the nuclear thickness function since it cannot be afforded to frequently perform 
time-consuming calculations. The parametrisation used here is (Fig. 110. 2j) 

^^-^'^ = I ; otherwise ^^^'^^ 

The radius of the nucleus Rau — 6.86 fm is taken as the zero of the polynomial. 
The parameters are 



ao = (2.193 ± 0.13) fm-2 

02 = (-0.0157 ± 0.0018) fm"^ 

a4 = (-0.000654 ± 0.000046) fm"*' 

The hard scattering events have to be distributed over the overlap region. Since 
the two nuclei are strongly Lorentz contracted their longitudinal extension is ne- 
glected so that all hard scattering points lie in the z = plane. The instant of 
overlap between the two nuclei also defines the time origin, i.e. all hard scatter- 
ings occur at t = 0. It is reasonable to assume that in the transverse plane the 
probability u{x,y) for hard interactions is largest in the region where the number 
of binary nucleon-nucleon collisions is highest. This implies that u{x,y) is coupled 
to the nuclear thickness function: 
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Glauber model calculation ^ 




Figure 10.2: Glauber calculation j^H] and parametrisation llU.41 for the thickness 
function of a gold nucleus 



u{x,y) oc Tau{x - b/2,y) ■ Tau{x + b/2,y) (10.5) 

Consequently the hard scatterings take place preferentially near the centre of the 
overlap region. 

There is a possibility to define a transverse area for the QGP which is smaller 
than the overlap region. Consequently not all hard scatterings take place inside 
the QGP volume and there is a finite propbability for these partons to escape 
undisturbed. 

10.1.2 Cronin Effect 

The Cronin effect is taken into account by a broadening of the intrinsic /c^-distribution 
of the partons (Ch. 15. 5j) . In this model the protons experience rescattering before 
the hard scattering. The additional transverse momentum built up by the protons 
translates into a bigger intrinsic transverse momentum of the partons: 

^fcx(^'^'^) = ^fc±,o + « ■ iNco\\{x,y,b) - 1) (10.6) 

where (Jkj^fl = 1 GeV is the width of the /c^-distribution used in proton- proton 
collisions and a is a constant. Here does not depend on Q, since parton showers 
are done explicitly in the simulation. This includes the evolution down to a fixed 
scale Qq and the intrinsic k± only has to cover the range below that. 

The number of collisions that a proton experiences before the hard scattering is 
estimated again with the help of the Glauber model. It depends on the point in 
the transverse plain where the proton is located (and where the hard scattering will 
happen) and the impact parameter. The maximum number of these rescatterings 
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is given by 

N^^{x,y,b) = aZ,TAu{x ± b/2,y) (10.7) 

with (T?fg[ ^ 42 mb at y/s = 200 GeV. The sign in the argument of Tau depends on 
the membership of the nucleon under consideration in one or the other nucleus. 

For simphcity it is assumed that the probabihty for a hard interaction in a binary 
nucleon-nucleon coUision is Phs and that there is no possibihty that the pair doesn't 
interact at alL Then the probabihty for the nucleon to make a hard interaction 
with the nth nucleon from the other nucleus is 

P„ = Phs(l-Phs)"-' (10.8) 

This is used in the model to choose for each nucleon the actual number of rescat- 
terings. Then the mean of the two values is taken because in Pythia the k±- 
distribution is always the same for both participants. For P^s the ratio of the cross 
section for hard scattering to the total cross section is used so that Phs ~ 0.04. 

The value of a can be determined from the dAu data. In that case only the 
nucleon from the deuteron has undergone rescattering and the one from the gold- 
nucleus is unchanged. This situation is equivalent to a symmetric increase with 
a/2 instead of a on both sides. For the comparison to data standard proton-proton 
collisions were simulated with the /c^-broadening as the only change. 

10.1.3 Model of the QGP: Energy Density Distribution, EOS 
and Evolution 

Since it is not a priori clear what the energy density distribution should be, there 
are two options available in the model: One is a constant, i.e. the energy density 
depends only on the proper time r, and in the other e is a function of the transverse 
coordinates x and y and r. In the latter case it is assumed that the energy density 
is proportional to the number of binary collisions, i.e. 

e{x,yAr = 1 fm/c) = Af ■ Tau(x - b/2,y) ■ Tau(x + b/2,y) (10.9) 

It is plausible that the energy deposition and thus the resulting energy density is 
larger in areas where more particles collide. The normalisation is chosen such that 
the mean energy density for 6 = and at r = Ifm/c is equal to the initial energy 
density. 

(Rau \ 

As a consequence the energy density is higher in the centre than near the periphery 
(Fig. 110. Hj) and on average decreases with increasing b. In this scenario it is natural 
to have the contour, where the energy density falls below the critical value, as 
boundary of the QGP and not the overlap region or some other artificial choice. 
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Figure 10.3: Energy density distribution in the transverse plain at r = Ifm/c for 
6 = 4fm 



The equation of state of an ideal ultra-relativistic gluon gas is used (Sec. 15. 2j) . 
The contribution from quarks is neglected, because the plasma is gluon-rich in the 
beginning due to the large gluon content in the proton at the typical x of interest at 
RHIC. Besides the gluons come into equilibrium very fast, because the gluon cross 
sections are large. Much more time is needed for the quarks. The crucial quantity 
is the gluon number density which is given by 



n = 4T3C(3) = ^^T=^ (10.11) 



Another useful relation is the one between the temperature and the energy density 
(Eq.EIl: 



For the evolution of the QGP a Bjorken-like scenario with only longitudinal 
expansion (Sec. I5.3|l is used. With the initial condition e(ro) = eo the energy 
density follows the evolution (Eq. I5.36j) 

<^) = ^o(^^J (10.13) 
while the temperature changes according to 

30 6„^'" 



This relation can be used to specify the lifetime of the QGP, defined as the point 
when the temperature drops below the critical temperature T^. 
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(10.15) 



The time evolution of the gluon density is correspondingly given by 

16-1.2 /30eo^^^'^ 



= — ^-l 77^) ■(-) (10.16) 



The plasma gluons have a mean thermal energy of 2.7 T which is in practice 
smeared out in the form of a Gaussian distribution with width ag. This motion is 
undirected. In addition there is the directed expansion motion in ±z direction with 
velocity z/t. Furthermore the gluons have a small effective mass rrig to account for 
interactions. 

The rapidity density of QGP gluons can be estimated easily if the thermal motion 
which is undirected and will cancel on average is neglected. Then there is a simple 
relation between the rapidity y and (here the spacial coordinates are denoted 
with Xi, X2 and X3 to avoid name conflicts with the rapidity): 




2 \P^ — 'P\\J 2 \p'^—p^xs/t 
which can be solved for X3 leading to 

X3 = Tsmh{y) (10.18) 
The rapidity density of plasma gluons at mid-rapidity is found to be 



diV„ 



dy 



y=0 



dxidx2 

7lRl^n{T)T 

993 



d^N dx3 



da:ida:2da:3 dy 



y=0 



(10.19) 



where it was assumed that the energy density is constant over the transverse plane. 
This is consistent with the numbers derived in |67] . 



10.1.4 Energy Loss 

The soft interactions (SI) of a traversing energetic parton with the QGP are treated 
as elastic scatterings from the plasma gluons. These are non-perturbative interac- 
tions and their characteristics can therefore not be derived from theory. So elastic 
scattering is the natural ansatz for such interactions between pointlike partons. It 
is assumed that the individual scatterings are independent. The parton interacts 
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with each gluon on its way out with a probabihty Pint- -Pint can be inherited from 
the SCI model where it was found to be roughly 0.5 (Ch. 

A suitable variable to describe the scatterings is the Mandelstam variable t, taken 
to be Gaussian distributed since the matrix elements calculated in the framework 
of pertubation theory are not applicable in this regime. The width at and mean nt 
can not be larger than a few hundred MeV for the process to be non-perturbative. 
The t-range is limited by the available energy: 

h = 4:{ml -E^)<t<0 = t2 (10.20) 

where nip and Ep are the parton's mass and energy in the cm. system. The 
scattering angle in the cm. system is given by 

In the lab system the parton is typically much more energetic than the plasma 
gluons so that it loses energy due to the boost and suffers only a small deflection. 

Gluons (as parton) should interact more strongly than quarks since they carry 
colour octet charge. In the model this is taken into account by giving the gluon a 
second chance if it didn't interact with a certain background gluon in the first try. 
Two interactions with the same gluon are excluded so that the effective interaction 
probability for a gluon becomes Pf^^ = 0.75. 

The parton is sensitive to a cylindrical volume with radius -Rcyi around its track, 
i.e. it has the possibility to interact with all gluons within this volume. The total 
number of encountered gluons is given by 

'^out g 

A', = >riJ?„/drn(r)=.fl?,,15_i£(^) (io.22) 

Tin 

where Tin and Tout are the proper times when the parton enters and leaves the 
QGP. In a scenario where the QGP fills nearly the entire overlap region Tin is 
simply the QGP formation time Tj, otherwise the geometry has to be considered as 
well. Similarly Tout is the point when the parton geometrically leaves the QGP but 
it is limited by the lifetime of the plasma: 



Tin > Ti (10.23) 
Tout < Tf (10.24) 

The interactions with the QGP are simulated explicitely in an iterative way: 
Starting from the entrance point the position of the next QGP constituent is es- 
timated and from a Gaussian around that point a choice is made for the actual 
position. Then the gluon is generated with appropriate thermal energy and expan- 
sion velocity and the soft scattering is simulated. The procedure is repeated going 
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out from the interaction point until Tout is overstepped or - if the energy density 
distribution is inhomogeneous - the local energy density falls below the critical 
value. 

It is not clear, how the proton remnants should be treated. Luckily, the ob- 
servables studied here are not sensitive to the treatment of the remnants, since 
the analysis is performed only in the central rapidity region and the remnants do 
not contribute to the high-p^ particle yields at mid-rapidity. Therefore, they are 
neglected and do not interact with the plasma at all. 

10.1.5 Hadronisation 

The interactions with the background also lead to changes in the string topology. 
The standard Lund string topology in shown in Figure 12.41 The interactions with 
the QGP change the strings and the partons presumably get connected to the 
plasma instead. It is also important to note that gluons radiated from a scattered 
parton in final state parton showers do not any more belong to the same string as 
the original parton. 

However, it is not clear how this system hadronises and what happens if the 
formation time of Ifm/c in the particle's rest frame lies within the QGP. The 
matter clearly calls for detailed investigations. For the study presented here a 
simple interim solution is employed in the SCI jet quenching model: The partons 
are hadronised with independent fragmentation (Sec. 12.3^ . This is sensible insofar 
as the parton will in any case have lost the connection to the proton remnant in the 
interactions with the background. The non-conservation of energy- momentum is 
not a problem in this context because the QGP acts as a kind of 'energy reservoir'. 

The pp reference and dAu runs are hadronised with Lund string fragmentation 
since there are no soft colour interactions in these cases. Strictly speaking parton 
systems that are created outside the QGP and do not interact at all would also 
have to be hadronised with the Lund model, however, this is not done for practical 
reasons. The effect would be very small anyway in central and mid-central collisions, 
because the transverse overlap area outside the QGP is small compared to the QGP 
area and scarcely populated due to the concentration of hard scatterings towards 
the centre. The situation might be different in peripheral collisions, depending on 
the size of the possible hadronic layer around the plasma. 

10.1.6 Model Parameters 

The parameters of the model are given in Table IIU.II together with their default 
values. The size of a hadronic layer can only be given when the energy density 
distribution is fiat, otherwise the size of the QGP is defined by the contour of the 
critical energy density. The values of parameters like the QGP formation time have 
been constrained in independent investigations so that they are not free, but can 
be varied in a certain range only. 
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The parameter for the Cronin effect is fixed with the help of dAu data. This is 
independent of the rest of the model, since in dAu collisions no QGP is formed and 
the part that describes the plasma and the interactions with it is irrelevant. The 
only difference to ordinary pp collisions is the /c^-broadening. 

The interaction probability is taken from the original SCI model and is not 
available as a free parameter in the jet quenching extension. 

The width of the t-distribution is taken to be constant but scaled down for 
parton-gluon pairs with small invariant mass. This is done to prevent the routine 
from becoming very inefficient when a lot of values have to be rejected because they 
lie outside the allowed range (see Eq. Il().2()j) . 

Table 10.1: Parameters of the SCI jet quenching model, their default values and 



constraints 
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1 r 1 j_ 1 

default value 
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U im 
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(only with homoge- 










neous energy density) 










QGP formation time 




0.2fm/c 


0.2. ..Ifm/c lEHl 




initial energy density 




5.5GeVfm-3 


4. ..6.5GeVfm-3 [HH] 




(to = Ifm/c) 










critical temperature 


T 


0.175 GeV 


0.17... 0.18 GeV 




width of gluon energy 




0.1 GeV 


? 




distr. 










gluon mass 


nig 


0.2 GeV 


? 


SI 


interaction probabil- 


Pint 


0.5 


~ 0.5 




ity 










radius of parton cylin- 




0.3 fm 


~ 0.3 fm 




der 










width of t distr. 




min(0.5GeVMti|) 


0.1...0.6GeV2 (?) 


Cronin 


increase of al per 


a 


0.25 GeV^ 


? 


effect 


scattering 









10.2 Results 

The event analysis concentrates on the following parts: 

• Ratio -RauAu between the p^-spectra for charged particles and neutral pions 
with |?7| < 0.35 and a pp reference 

• Centrality dependence of -RauAu 

• Azimuthal correlation of charged particles with |?7| < 0.7: Particles with 
4 < < 6 GeV are defined as trigger particles, for each trigger particle the 
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azimuthal separation A0 from all other particles with 2 < p± < p_L(trig) 
in the event is calculated. In the analysis considered here several trigger 
particles in one event are also possible, but this happens only in ~2% of all 
events containing trigger particles. 

• Change of the azimuthal correlation with respect to the reaction plane orien- 
tation: trigger particles with (p^"^ < ir/A and (p^"^ > Sir/A are defined to be 
in-plane and those with n/A < 0*'''^ < Sir/A are out-of-plane 

The pp reference was generated by running standard pp in the same event generator. 
The Cronin effect is the only nuclear effect that is explicitly taken into account. 
One difference between the measured and the model nuclear modification factor is 
the effect of different A-scaling. In the data this leads to a decrease at small p± even 
in the absence of all nuclear effects which is not present in the model. However, 
the absence of this scaling-effect in the simulation does not influence the results 
because the model is not applicable to the low-p± regime anyway. 

The azimuthal correlation should be the same as the one measured by the STAR 
collaboration after their subtraction of the elliptic flow and pedestal contribution 
(collective flow is not included in the model). 

The particle production mechanism has not been modified which means that the 
relative yields of hadron species are unaltered. In particular there is no excess of 
(anti)protons over pions as observed in experiment (Ch. 16. 2j) . The model results 
should thus in the first hand be compared to the 7r° data when possible. In the 
model the effects are the same for all particle species and the ratio is calculated 
for the sum of the spectra for neutral pions and charged particles in order to get a 
better statistics. 

There are several reasons why the model cannot be considered trustworthy for 
P± ^ 3 GeV. On the parton level only hard processes with p_L > 3 GeV are gener- 
ated, all softer (perturbative and non-pertubative) interactions as well as multiple 
interactions that affect the low-p_L region are not taken into account. In a strict 
sense the p±-cut on the matrix element influences the p_L-spectra of hadrons up to 
4-5 GeV: Because of the intrinsic k± too few particles are produced in this regime. 
This effect will cancel only partly in the -RauAu since the effect of /c^-broadening is 
not present in the pp reference. 

The p±-cut has practically no influence on the azimuthal correlation since a 
trigger particle with p± > A GeV is required. Then the hard scattered parton also 
had a. p± > A GeV (larger than the hadron's p±) and thus stems from a region that 
is hardly affected by the p±-cut. In the framework of the SCI jet quenching model 
all final state particles originate from the fragmentation of the high-p^ parton, this 
process is entirely independent from the hard scattering and cannot be affected by 
the cut-off. 

Another aspect that is not included in the model is the hadronisation of the 
plasma. This is also expected to contribute to the low-energy particle yields. 
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All results shown in this section were obtained with the default values of the 
parameters given in Table ITim "SCI jet quenching model I (II)" refers to the 
ansatz with constant (inhomogeneous) energy density distribution. 

The statistical uncertainty on the model result for the nuclear modification factor 
increases with p±: It is < 5% below 6GeV, ~10% below 8GeV and > 20% above 
8 GeV. On the azimuthal correlations it is smaller, namely of the order of a few 
percent. 



10.2.1 Cronin Effect 

The parameter a for the Cronin effect was fixed by a comparison to the dAu 
data. This was done by running the simulation with the fc^-broadening as only 
modification to standard pp. In Pythia the k_i_ is always the same for both colliding 
partons, but in dAu only one nucleon (that from the deuteron) has experienced 
rescattering. This situation is mathematically equivalent to a symmetric increase 
of the variance with a/2 instead of a (the Gaussian distributions of the two protons 
are convoluted leading to a Gaussian whose variance is given by the sum of the 
variances of the initial distributions). The effect of the second nucleon in the 
deuteron has been neglected. In fact is has been shown experimentally in |1U] that 
RdAn and -RpAu are very similar. 
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Figure 10.4: Model results for dAu (charged particles and neutral pions at |?7| < 
0.35) and PHENIX vr^ data [lU] (errorbars indicate statistical and 
boxes systematic errors) 

The best results were obtained with a = 0.25 GeV^ and are shown in Figure ITO. 41 
together with the experimental data. The agreement is reasonably good for 3 < 
p_L < 9 GeV, i.e. where the p± is large enough for the model to be reliable and not 
too large so that the statistical error is moderate. 
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There is, however, a problem with the value of a = 0.25 GeV^: One might want to 
interpret it as the square of the mean momentum transfer per scattering because of 
the way the model is constructed fCh. 15.51 and llO.l"^ . But the momentum transfer 
should not exceed 0.2 GeV if the proton is supposed to interact as a whole. This 
number is motivated by the uncertainty principle: The size of a proton is roughly 
1 fm and a resolution of 1 fm^^ corresponds to a momentum transfer of 0.2 GeV. a 
being considerably larger than 0.04 GeV^ can be seen as evidence for other effects 
that build up additional transverse momentum. However, also a highly excited or 
dissociated proton interacts so that a larger value of a may be reasonable. The 
model for the Cronin effect used here is then an effective description that should 
not be taken too literally. 

A variation of a by ±0.1 GeV^ leads to only small changes with the result still 
being consistent with the data. 

Figure ITO . 51 shows that the broadening of the intrinsic also has a small influence 
on the 2-particle azimuthal correlation. The fraction of events containing trigger 
particles is slightly increased from 3.16 ■ 10~^ in pp to 3.75 ■ 10~^. 
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Figure 10.5: 2-particle azimuthal correlation in the SCI jet quenching model for 
pp and dAu collisions 



10.2.2 Nuclear Modification Factor 

The model results for the nuclear modification factor are shown in Figure 110.61 
together with the experimental data for five different centrality classes. In the 
most central and most peripheral class the effect of the /c^-broadening without 
scatterings is shown in addition. It naturally increases with centrality because 
the nuclear overlap is much larger in central collisions. Thus the Cronin effect 
counteracts the partonic energy loss which is larger than it seems at first sight. 
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Figure 10.6: Model results for the nuclear modification factor with constant en- 
ergy density for different centrality classes and PHENIX data j3T], 
for 0-10% and 80-92% centrality also the effect of only fc^-broadening 
without partonic energy loss is shown 



66 



The approximately flat shape of -RauAu for p± > 3 GeV is well reproduced by 
the model. There is a slight increase in the range between 3 and 4-5 GeV that is 
caused by the p±-cut on the matrix element. Because the intrinsic k± has to be 
added to the momentum transfer in the hard scattering, there are too few particles 
produced also above the threshold. The effect is more pronounced in nucleus- 
nucleus collisions due to the Cronin effect and thereore it doesn't cancel completely 
in the ratio. As expected it increases with centrality and also reaches out to higher 
Pj_ in more central events. 

In the most peripheral class the model result is in agreement with the data but 
is somewhat low. The situation is slightly better in 60-70% but already in 40-50% 
the model result falls behind the data which decrease faster. The trend advances in 
the two central classes so that in the 0-10% bin the model gives roughly 50% of the 
effect observed in the data. The centrality dependence of the model result is rather 
weak so that it is practically impossible to identify the functional dependence. 

With the non-uniform energy density distribution the result for 0-10% centrality 
is nearly identical with the one obtained with a constant energy density. But the 
increase when going to more peripheral collisions is stronger because the mean en- 
ergy density decreases. It thus agrees better with the data than the other approach. 
This is also in reflected in Figure ITO . 81 where it becomes clear that the model follows 
the data for peripheral events, but then it levels off and stays approximately con- 
stant. In fact, the model depends quadratically on the centrality, there is no linear 
contribution. As will be discussed in Chapter ^2 this is related to the dependence 
of the energy loss per interaction on the parton energy. Figure 110.91 shows the 
energy loss in the interaction with the first plasma gluon (interaction probability 
is 1 in this case) for light quarks as a function of the initial quark energy. Even 
for the highest energies the energy loss rises slightly (the averaged values after 10 
scatterings are below those after only one interaction). On the other hand, quarks 
with very small energy gain energy because the thermal energy of the environment 
is higher than the quark energy. 
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Figure 10.7: Model results for the nuclear modification factor with non-uniform 
energy density for different centrality classes and PHENIX data jM] 
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Figure 10.8: Comparison of the centrality dependence of the model with inhomo- 
geneous energy density profile to PHENIX data j3J 
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Figure 10.9: Dependence of the energy loss per scattering on the initial energy 
of the quark for u and d quarks, scattering centre is the first gluon 
(" 10 scatterings" denotes 10 interactions with the same gluon) 
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10.2.3 Azimuthal Correlation 
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Figure 10.10: 2-particle azimuthal correlation from the model with constant (up- 
per figure) and non-uniform energy density profile (lower figure) in 
peripheral (80-92%) and central collisions (0-10%) 

The 2-particle azimuthal correlations for the two versions of the model are shown 
in Figure 110.101 There is a suppression of the away-side jet but no disappearance 
as in the data. Again the ansatz with an inhomogeneous energy density profile 
leads to a bigger effect. But while the Cronin effect counteracted the energy loss in 
the case of the nuclear modification factor, it now contributes to the jet quenching 
fFig. llO.TTj) . It leads to an away-side peak that is lower in central than in peripheral 
events, although it does not help on average fFig. 110. 

The near-side jet remains unchanged, but as discussed aheady in Section 16.2.21 
this does not necessarily mean that the respective parton has not suffered energy 
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Figure 10.11: Azimuthal correlation in AuAu collisions with only Cronin effect 
(no scattering) for central and peripheral events 



loss. Since the partons are assumed to hadronise outside the QGP and ordinary 
vacuum fragmentation functions are used, the shape and the integral of the near- 
side peak are determined only by the momentum of the partons. The energy loss 
in the QGP shifts the parton p^-spectrum to smaller momenta so that the overall 
yield of partons, that have enough energy to produce a trigger particle, is smaller. 
But since the 2-particle azimuthal correlation is normalized to the number of trigger 
particles and contains no information about the rate at which these particles are 
produced, this effect is invisible. However, there is a certain sensitivity to the shape 
of the parton p^-spectrum because a more energetic parton will on average produce 
more associated particles. A flat spectrum will thus lead to a near-side peak with 
a larger area than a steeply falling spectrum because the relative contribution from 
partons with very high momenta is smaller in the latter case. 

It is thus oversimplified to conclude that the partons producing trigger particles 
have not lost a substantial amount of energy and therefore stem form the periphery 
of the QGP. In fact, it can easily be shown that this is not the case. Figure ITO. 121 
shows the results from a simulation with the SCI jet quenching model II where all 
hard scattering points were placed in a 0.5 fm thick shell at the edge of the overlap 
region (b = 0). At first sight it is surprising that the jet quenching is very weak in 
this configuration but it can be explained as a mainly geometrical effect: If a parton 
is emitted from the surface with random direction in the transverse plane (this is 
the case for hard scattering) only a range of ~ 7r/2 in angle will lead to a sizable 
path through the plasma. Furthermore, the constituent density is much smaller 
near the edge than in the centre. But while a parton directed towards the centre of 
the QGP moves inward, the number density drops rapidly so that the parton never 
reaches a region of high density. The lifetime of the QGP is, with the parameters 
chosen here, rj ~ 5fm/c so that it hadronises even before the parton has reached 
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the centre. The sum of these effects leads to the very weak jet quenching observed 
in Figure ITO. 121 The jet quenching increases with the thickness of the shell until for 
a thickness that equals the nuclear radius the original result for central collisions 
is recovered. The result is very similar for a constant energy density profile. This 
means that the observed hadrons are the fragmentation products of partons that 
come from practically everywhere in the overlap region. 
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Figure 10.12: Results of the simulation with the SCI jet quenching model II but 
now all hard scattering points lie in a narrow shell at the edge, 
centrality is 0% 



To get a substantial suppression of the away-side jet one parton has to lose much 
more energy than the other one. Since the crucial quantity for energy loss is the 
number density of constituent gluons, a situation where one parton encounters 
much more gluos is needed. There are several reasons why this is hard to achieve 
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(and thus the effect on the azimuthal correlation is so small). Firstly the hard 
scattering points are concentrated towards the centre of the QGP and consequently 
the difference between the path lengths in the plasma of the two scattered partons 
is typically small. What makes the situation even worse is that the gluon density 
falls like r~^. In fact, the energy loss per scattering increases as the temperature 
falls but the decrease in number density prevails (the consituent gluon energy is 
proportional to t~^/^ only). This means that most of the energy loss happens at 
early times and small differences in path length are further reduced because the 
plasma is already quite dilute when the first parton leaves the overlap region. It is 
also seen in Figure IIO.IH^ which shows the specific energy loss of light quarks as a 
function of the path length in the QGP. Quarks with low momentum gain energy 
because they are below the thermal energy in the beginning. They are also more 
sensitive to the plasma gluon energy than more energetic quarks, i.e. they follow 
the evolution of the QGP more closely. 
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Figure 10.13: Specific energy loss of monoergetic light quarks emitted from the 
centre of the QGP with ^ = -7r/4, energy density is homogeneous 

Furthermore, the QGP lifetime limits the possible path length differences in 
central collisions. Even for parton pairs that happen to be produced at the surface 
the path length difference cannot exceed the lifetime of ~ 5 fm. The model with 
the non-uniform energy density profile works better because it helps to amplify 
differences. If the hard scattering occured not directly in the centre the parton with 
the longer path traverses the denser central region and experiences more interactions 
than the other one which only moves through the more dilute outer parts of the 
QGP. 

The asymmetry between the two hard scattered partons is enhanced by the ma- 
trix element since ~50% of all parton level processes are q + g ^ q + g scatterings. 
The gluons interact more strongly and lose more energy than the quarks. 
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Figure 10.14: Correlation of tlie 2-particIe azimutlial correlation with the reaction 
plane from the model with constant (upper figure) and non-uniform 
energy density profile lower figure) in mid-central collisions (20- 



In Figure ITO. 141 the correlation with the reaction plane is shown for mid-central 
events. In both cases there is no difference between the particles emitted in-plane 
and out-of-plane. The reason is that this observable suffers even more from the 
problems with the similar path lengths. The mean radius of the overlap region is 
3.0 fm in-plane and 4.3 fm out-of-plane for 40% centrality. The difference is at most 
1.3 fm in the whole range from 20 to 60% centrality so that the differences in path 
length of the two partons are small. Due to the other difficulties discussed before, 
it is not surprising that there is no visible effect on the azimuthal correlation. 



75 



10.3 Discussion 



The SCI jet quenching model yields good results for the nuclear modification factor. 
The shape is similar to the data, namely approximately fiat for p± > 3 GeV and 
the magnitude can reach 50% of the observed effect in central collisions. The rest is 
presumably due to gluon bremsstrahlung which should naturally come along with 
the scatterings. The centrality dependence is, however, not linear as in the data 
but quadratic. It is not clear how this will change when gluon bremsstrahlung is 
included. 

There is also a clear suppression of the away-side jet but the effect is much 
smaller than observed in data. The reasons for this are mostly understood: The 
hard scattering points are concentrated towards the centre of the QGP, the gluon 
density in the plasma falls rapidly and the lifetime limits the available path length 
differences between hard scattered partons. These are not details of this model but 
very general features. 
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Figure 10.15: 2-particle azimuthal correlations in pp collisions for different frag- 
mentation schemes 



Another aspect is that the azimuthal correlation depends on the details of the 
fragmentation scheme. This is illustrated in Figure (TO . 1 51 which shows a comparison 
of Lund string fragmentation to independent fragmentation in pp simulations. The 
fragmentation function is the same in both cases (namely the Lund symmetric 
function). The independent fragmentation scheme produces somewhat less particles 
at midrapidity than the string fragmentation, but this is a small effect (the ratio of 
the p±- spectra is ~ 0.95). For the 2-particle azimuthal correlation, however, this 
makes a big difference. Modifications of the fragmentation scheme and function 
are a point that has not been addressed in detail in this model but it is propbable 
that the presence of a QGP will also infiuence the fragmentation. This may lead to 
substantial changes in the azimuthal correlations, a softer fragmentation function 
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for instance will enhance the suppression of the away-side jet. 

Still it can be concluded that it is indeed possible to obtain a considerable contri- 
bution to the observed jet quenching from energy loss due to soft scatterings. The 
scenario with an energy density that is higher near the centre that at the edge of 
the QGP yields better results. 

The energy loss per scattering is governed by the width of the t distribution 
at and, to a smaller extent, the gluon mass rrig. A more massive gluon leads to 
a smaller energy loss. These parameters affect the overall energy loss and the 
centrality dependence. 

The QGP formation time, initial energy density and critical temperature de- 
termine the number of scatterings but don't affect the centrality dependence. A 
shorter formation time for example will add a certain number of encountered gluons 
to each parton irrespective of the centrality. Because of the high energy density in 
the early stage the overall energy loss depends quite strongly on the formation time 
and moderately on the initial energy density but hardly on the critical temperature. 

The Cronin parameter is fixed independently with the help of the dAu data 
but there is still some freedom. A stronger fc^-broadening reduces the energy loss 
such that the centrality dependence of -RauAu becomes flatter but it increases the 
suppression of the away-side jet in central collisions. 

The effects of a hadronic layer around the plasma phase have not been studied 
here, but their overall characteristics can be qualitatively described. It would de- 
crease the energy loss, but not very much as long as the hadronic layer is thin. This 
may in turn alter the influence of the centrality, depending on how the thickness 
changes with it. A more important effect could enter in connection with modi- 
fied fragmentation. Partons that originate from scatterings in the hadronic layer 
may (depending on their direction of motion) escape without any interactions and 
hadronise as in vacuum. This has a potentially sizable impact on the azimuthal 
correlation. 
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11 Exploring General Properties of 
Partonic Energy Loss 



11.1 Motivation 

It has been shown in Section 15.61 that the dependence of the energy loss AE on 
the parton energy E and the path length L (in a static medium) is sensitive to 
the physics of the energy loss mechanism. Matter induced gluon radiation can for 
instance lead to a typical L^-dependence, if there are coherence effects. The energy 
loss by incoherent scatterings on the other hand is proportional to L. It is thus 
obvious to ask if it is possible, with the presently available data on jet quenching, 
to distinguish between different scenarios. 

Now it is fortunate that the SCI jet quenching model is constructed in a modular 
fashion. It is thus possible to keep the model for the QGP but replace the energy 
loss by something less specific, so that different dependencies can be tested and 
compared to data. This can of course not replace a physical model, but it can easily 
be used to investigate if the data favour a certain type of energy loss mechanism. 

11.2 Energy Loss in the Toymodel 

The only part of the model that is different from the SCI jet quenching model is 
the energy loss. It is not simulated explicitely in the framework of the Toymodel, 
instead a certain amount of energy is substracted for each parton and deflection 
is not taken into account. This very general ansatz makes it possible to explore 
different dependencies of the overall energy loss. Since the path length is not a 
meaningful quantity in an expanding QGP, the number of plasma gluons Ng, that 
a parton encounters on its way, is used as a measure for the "amount" of plasma 
seen by a parton (this is sensible since also for gluon bremsstrahlung the number 
of scatterings is the important quantity). The dependence of the energy loss on 
Ng and on the parton energy can hence be varied according to the different cases 
listed in Table ITTTTl The constant of proportionality is a free parameter that has to 
be adjusted. Again gluons lose more energy than quarks, in the Toymodel this is 
achieved by doubling the parameter. If the energy of a parton drops below a certain 
value /ith + o"th, it is assigned a new momentum vector with energy picked from a 
Gaussian with mean /ith and width ath, the direction of motion is chosen randomly. 
This is of course a very rough way to model the thermalisation of partons, but since 
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they are anyway far below the trustworthy region in p± of the model the details of 
the process are not important. 

Table 11.1: Versions of energy loss in the Toymodel, energies are in units of 1 GeV 



Name 


Dependencies 


Constant 




path length 


energy 


energy density 




Toymodel laa 




E 


uniform 


0.065 GeV 


Toymodel lb a 


^1 


E 


uniform 


0.016 GeV 


Toymodel II a a 




E 


non-uniform 


0.070 GeV 


Toymodel lib a 




E 


non-uniform 


0.018 GeV 


Toymodel II a /3 




^/e 


non-uniform 


0.300 GeV 


Toymodel II b /3 




4e 


non-uniform 


0.020 GeV 


Toymodel II a 7 




const. 


non- uniform 


0.425 GeV 


Toymodel lib 7 




const. 


non-uniform 


0.200 GeV 



The other parameters that describe the QGP are the same as in the SCI jet 
quenching model. A complete list with the default values is given in Table 111.21 
The fragmentation scheme is again independent fragmentation. 



Table 11.2: Parameters of the Toymodel, their default values and constraints 





Parameter 


default value 


sensible range 


QGP 


size of hadronic layer 
(only with homoge- 
neous energy density) 


d 


Ofm 


0. ..3fm (?) 




QGP formation time 




0.6 fm/c 


0.2. ..1 fm/c ^ 




initial energy density 




5.5GeVfm-3 


4. ..6.5GeVfm-3 ^ 




(ro = Ifm/c) 










critical temperature 


T 

^ c 


0.17GeV 


0.17... 0.18 GeV 




radius of parton cylin- 


-Rcyl 


0.3 fm 


~ 0.3 fm 




der 








Cronin 


increase of a^^ per 


a 


0.25 GeV2 


? 


effect 


scattering 








energy 


mean thermal energy 




0.6 GeV 




loss 


width of th. energy 
distr. 

constant of propor- 
tionality 




0.1 GeV 





The constant is adjusted in one centrality bin (either the most central or the 
most peripheral) and then the other classes are generated with the same value. 
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11.3 Results 



In this section only a fraction of the results is shown, the complete set can be found 
in Appendix El The shape of the nuclear modification factor is very similar in all 
versions: It is nearly flat with a tendency to decrease slightly towards high pj_, only 
when AE is independent of E it rises. The centrality dependence and the azimuthal 
correlation are more interesting, because they are sensitive to the differences in the 
model assumed for the energy loss. 
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Figure 11.1: Centrality dependence of -RauAu for the Toymodel laa with AE oc 
NgE and Toymodel Iba with AE oc N^E, both with uniform energy 
density profile 

Fi gure fTl . 1 1 shows the centrality dependence of -RauAu for the configurations with 
a homogeneous energy density profile Toymodel laa and Iba. The shape is in 
both cases similar, although the linear version (I a a) seems to be somewhat flatter. 
The course of laa is mildly curved while Iba is more straight. Both model re- 
sults are significantly flatter than the data. The absolut magnitude of the nuclear 
modification factor is not meaningful, since it is governed by a free parameter. 

The effect on the 2-particle correlation is in both cases small, but stronger in Iba 
(Fig. So altogether the configuration with AE oc Ng appears to be slightly 

closer to the experimental data. 

In the configuration with a non-uniform energy density the mean energy density 
increases with centrality. Thus the expectation is that the shape of the centrality 
dependence of -RauAu is steeper than with a constant energy density. This is in fact 
observed (Fig. I11.3jl . Although both calculation results are slightly curved, they 
now give a reasonably good description of the data. The bending is again more 
distinct in the linear ansatz. 
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Figure 11.2: 2-particle azimuthal correlation for tlie Toymodel laa with AE oc 
NgE and Toymodel Iba with AE oc N^E, both with uniform energy 
density profile, the 80-92% data are the same in both cases 
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Figure 11.3: Centrality dependence of -RauAu for the Toymodel II aa with AE oc 
NgE and Toymodel lib a with AE oc N^E, both with non-uniform 
energy density profile 



The suppression of the away-side jet has also improved in both compared 
to Toymodel I fFig. lll.4|) . This matches with the observation made with the SCI jet- 
quenching model where also the inhomogeneous energy density distribution yields 
better results. 

When AE oc the centrality dependence is parabolic (Fig. I11.5j] without a 
linear component and hence qualitatively different from the data. The two scenarios 
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Figure 11.4: 2-particle azimuthal correlation for the Toymodel II aa with AE oc 
NgE and Toymodel lib a with AE oc N^E, both with non-uniform 
energy density profile, the 80-92% data are the same in both cases 

1.2 I 1 1 1 



Toymodel llap 
Toymodel llbp 



> 
C5 

CO 

A 



0.6 



0.4 

V 



0.2 



I ->t r- 



20 40 60 80 100 

centrality [%] 

Figure 11.5: Centrality dependence of -RauAu ^ot the Toymodel llaj3 with AE oc 
Ng\fE and Toymodel IIb/3 with AE oc N'^\fE, both with non- 
uniform energy density profile 



Toymodel II a /? and Wh (3 are again very similar, but this time the curvature is 
somewhat stronger in the quadratic setup IIb/3. 

In the 2-particle correlation on the other hand the situation is as before: The 
suppression obtained with IIb/5 is clearly stronger than with IIa/3. The effect is 
even bigger than before and has reached a reduction by more than 50% fFig. lll.Bj) . 
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Figure 11.6: 2-particle azimuthal correlation for tlie Toymodel II a/? with AE oc 
Ngy/E and Toymodel IIb/3 with oc N^y/E, both with non- 
uniform energy density profile, the 80-92% data are the same in 
both cases 
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Figure 11.7: Centrality dependence of -RauAu for the Toymodel II a7 with /S.E oc 
Ng and Toymodel lib 7 with /S.E oc A'^^, both with non-uniform 
energy density profile 

Surprisingly the results of the ansatz with AE independent of E are very similar 
to those obtained with AE oc E. The centrality dependence is now linear in 
both cases and again the II a 7 is slightly flatter than lib 7. While the centrality 
dependence of -RauAu agrees quite well with the data the azimuthal correlation is 
not as convincing. There is a suppression which is again stronger for the quadratic 
version lib 7, but this is still far from a disappearance. It is also slightly worse 
than the results obtained with AE oc E. 
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Figure 11.8: 2-particle azimuthal correlation for tlie Toymodel II a7 with AE oc 
Ng and Toymodel lib 7 with AE oc Ng, both with non- uniform 
energy density profile, the 80-92% data are the same in both cases 



11.4 Summary and Discussion 

In summary, the differences in -RauAu and its centrality dependence between the 
configurations with AE oc Ng and AE oc A*"^ are small. AE oc Ng has the tendency 
to look slightly better and also the suppression of the away-side jet is stronger. 
However, these differences are not big enough to come to a clear cut conclusion. 
The Ng dependence seems to be favoured, but the Ng dependence is not ruled out. 

The situation becomes partly clearer when it comes to the energy dependence. 
Here the AE oc scenario leads to a centrality dependence of -RauAu that has a 
clear parabolic shape and is thus qualitatively different from the linear behaviour 
observed in the data. The suppression of the away-side jet is stronger than in 
the other versions, but since this observable is not very sensitive to the energy 
dependence this is not a strong argument. Besides, the azimuthal correlation is 
obviously quite sensitive to details of the fragmentation which is not well under 
control in this model. So it is justified to conclude that AE oc ^/E is clearly 
disfavoured. 

But then it becomes more difficult since it is practically impossible to distinguish 
between AE oc E and AE independent of E. Perhaps the centrality dependence 
is slightly better in the latter case, but for the azimuthal correlation it is just the 
other way around. There are minimal differences in the shape of -RauAu^ It has 
a tendency to fall mildly when AE oc E and rise when AE does not depend on 
E, but the data don't favour one or the other alternative (Fig. I11.9|l . With the 
presently available data it is therefore not possible to differentiate between the two 
options and also with better statistics it is not trivial since the differences in -RauAu 
are also at higher p± not big. 
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There is a difference between the two scenarios when going to higher p± with the 
2-particle azimuthal correlation. This is shown with an example for the AE oc A'^ 
configurations in Figure 111.101 and 111.111 The suppression of the away-side jet 
is clearly stronger in Toymodel lib a than in lib 7. It is also observed that the 
height of the away-side peak relative to the near-side peak increases when going 
to higher p_i_, although the nuclear modification factor stays more or less fiat. The 
increase is thus not due to a "punch-through", i.e. particles that have so much 
energy that they lose only a very small fraction and emerge nearly undisturbed 
from the plasma. In fact, it is connected to the observation that the jet quenching 
becomes inefficient when the hard scattering occurs near the surface of the QGP 
made in Section [10.2. HI When the trigger condition gets more restrictive only jets 
from partons that suffered hardly any energy loss are selected. This means that 
the hard scattering points of the selected events move more and more towards the 
surface and the quenching becomes more inefficient. Consequently the away-side 
jet reappears even though the nuclear modification factor stays constant. 
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Figure 11.9: -RauAu in Toymodel lib a with AE oc NgE and Toymodel lib 7 with 
AE oc 



The results of this section show that the gluon bremsstrahlung scenario in the 
high energy limit (Eq. I5.58|) is within the bounds of possibility while the result 
obtained for L > (Eq. seems to be ruled out, at least when it is not 

combined with some other mechanism. 

The SCI jet quenching model also has an energy dependence (Fig. 110. 9^ which is 
not consistent with the data. But here it is clear that there is an additional effect 
missing since it can account only for 50% of the observed jet quenching. And it 
cannot be expected that the SCI jet quenching model describes the whole effect. So 
there is still freedom to resolve the present difficulties when the model is combined 
with induced gluon radiation (or possibly even more effects) to get a full description 
of the energy loss. 
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Figure 11.10: 2-particle azimuthal correlation for tlie Toymodel lib a witli AE oc 
N^E and Toymodel lib 7 with AE oc for 0-10% centrality, 
trigger particles satisfy 6 GeV < p± (trig) < 8 GeV and associated 
particles have 4 GeV < p_L(assoc) < p_L(trig) 
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Figure 11.11: 2-particle azimuthal correlation for the Toymodel llba with AE oc 
N^E and Toymodel lib 7 with AE oc for 0-10% centrality, 
trigger particles satisfy 8 GeV < j>_L(trig) < 10 GeV and associated 
particles have 6 GeV < p± (assoc) < p± (trig 



The inhomogeneous energy density distribution works better than the homoge- 
neous. This agrees with the findings of the SCI jet quenching model and there is 
no reason why the energy density should be homogeneous. Concerning a possible 
hadronic layer the same considerations as for the SCI jet quenching model apply, 
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namely that it can have an impact on the influence of the centrality depending on 
how the thickness of the layer changes with centrality. 

The fact that the centrality dependence does not level off and become flat at a 
certain level illustrates that the available path lengths are limited predominantly by 
the geometry and not by the QGP lifetime. Otherwise, at some point an increase 
of the size of the overlap region (i.e. going to more central events) would not result 
in longer path lengths and the nuclear modification factor would be unchanged. 

Since the Toymodel makes use of the same model for the geometry (e.g. distribu- 
tion of hard scatterings in the overlap region) , the QGP etc. as the SCI jet quenching 
model, it is expected to suffer from the same problems with the 2-particle azimuthal 
correlations. However, because they are mainly due to very general properties of 
the QGP one would think that the energy loss mechanism has to compensate the 
difficulties. The fact that even those versions of the Toymodel, which yield a good 
description of the nuclear modification factor, fail to describe the disappearance of 
the away-side jet strengthens the case for the hypothesis that the azimuthal cor- 
relation is fairly sensitive to the fragmentation scheme and function. However, it 
seems that it is also influenced by the path length dependence of the energy loss 
mechanism. 

The nuclear modification factor on the other hand was found to depend only very 
weakly on the exact energy and path length correlation with the energy loss and 
is also insensitive to the energy density distribution. The centrality dependence of 
-RauAu probes mainly the energy dependence. 

Now one can also try to answer the question of how the linear centrality de- 
pendence comes about. Strangely enough, it seems that a lot of different effects 
conspire to yield a very simple form in the end. It has been shown that apart from 
the geometry (and the QGP lifetime) also the Cronin effect, the energy density dis- 
tribution, the energy dependence of the energy loss and to some extend the impact 
of the path length on the energy loss influence the behaviour with the centrality. 
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12 Conclusions and Outlook 



Given the wealth of different observables that point towards QGP formation at 
RHIC, it seems that indeed the long-sought quark-gluon plasma has been created. 
Among them, the jet quenching plays an important role. For SPS the question 
whether a QGP is formed or not is not so easy to answer. But as has been shown in 
this thesis, the shapes of the p^-spectra are not very different from those measured 
at RHIC. Since this applies also for other observables, it is very well conceivable 
that the phase boundary has already been crossed in central collisions at SPS. 
However, from SPS data it is very difficult to investigate the properties of the new 
state of matter. This is to some extend possible at RHIC, but still a very difficult 
task that requires a profound understanding in particular of the high-p^ physics 
and the jet quenching. 

It was shown that soft rescattering of a hard scattered parton in a QGP can 
contribute significantly to the overall energy loss. With the model introduced in 
the first part of this thesis, roughly 50% of the observed effect in central collisions 
was reproduced. The magnitude depends on the momentum transfer per scattering 
which is a parameter, but it is clear that losses induced by scattering may not 
simply be neglected. 

It is not a problem that the model cannot account for the whole effect, since it can- 
not be considered complete in the sense that all processes contributing to the energy 
loss are included. The most obvious example is medium induced gluon radiation. 
It seems natural that scattering of the parton should induce gluon bremsstrahlung 
resulting in further energy loss. Additional processes will of course also change the 
centrality dependence so that at this stage it should not be overemphasised that 
the centrality dependence found in the SCI jet quenching model does not match 
the data. 

Nevertheless, the model is in two respects not very well motivated. First, it could 
be argued that the assumption of independent scatterings is not justified, at least at 
the early times when the density in the QGP is high, since soft interactions happen 
on relatively long time scales. Then interferences might become important and it is 
not clear how they should be treated in a Monte Carlo generator. The second point 
is that since the momentum transfer is small it is impossible to resolve individual 
gluons. This would mean that the plasma has to be treated as a colour charged 
gluonic field with which the partons would nevertheless interact. But the parton- 
field interaction may not be very different from the parton-gluon interaction. Since 
there is poor theoretical guidance for soft interactions, the parton-gluon scattering 
can be used as a model. 
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It should also be noted that the (single) gluon exchange model should not be 
taken too literally. The Gaussian momentum transfer should rather be interpreted 
as due to a superposition of many-gluon exchanges. The exact details of this process 
are not relevant for this study. 

The main conclusion is unaffected by these problems that are unescapable when 
dealing with soft physics: Soft scattering can result in a sizable energy loss. 

The lesson to learn from the investigations of the overall features of partonic 
energy loss is that it is difficult to distinguish between different scenarios given the 
presently available data. AE oc seems unprobable because it yields the wrong 
centrality dependence of -RauAu- But AE oc E and AE oc const are practically 
undistinguishable, this will maybe become possible with better data at higher p^. 
AE oc Ng and AE oc Ng are also similar, AE oc A'^ looks slightly better. This 
may be interpreted as a hint for coherent gluon radiation but it is not a strong 
indication. In summary, it has not been possible to get to a clear cut conclusion 
for or against coherent gluon bremsstrahlung. 

The 2-particle azimuthal correlation seems to depend on details of the fragmen- 
tation procedure, an issue that has not been addressed in great detail in this study. 
This makes it very difficult to judge the results of the models but also to under- 
stand the data when attention is paid to the details. Beyond that, the azimuthal 
correlation shows a certain sensitivity to the path length dependence of the energy 
loss. In contrast the nuclear modification factor is insensitive to many details and 
the centrality dependence of -RauAu is ruled mainly by the energy dependence. 

For the future it is desirable to gain a better understanding of the interplay of 
scattering and gluon radiation. Unfortunately, including radiation in the SCI jet 
quenching model is not free of difficulties. While incoherent radiation is rather 
straightforward the coherence effects are extremely difficult to handle in a Monte 
Carlo framework. 

Another unsolved problem are the medium modifications to hadronisation. This 
is also highly non-trivial, mainly due to lack of information. Moreover, it is difficult 
to disentangle the effect of modified fragmentation from other medium modifica- 
tions. Hopefully, one can get some guidance from the observation of modified 
fragmentation in cold nuclear matter (e.g. 

If one wants to fully exploit jet tomography, i.e. use the jet quenching to de- 
termine the properties of the QGP, one first has to come to a consistent picture 
and get the partonic energy loss theoretically better under control. This means 
that further effort has to be made in order to gain a better understanding of the 
underlying physical principles. 
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Figure A.l: Fit parameters for the pion spectra 
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Figure A. 2: Enlarged low p± region of the above pion results 
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Figure A. 3: Fit parameters for the feed-down corrected proton spectra 
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Figure A. 4: Fit parameters for tlie feed-down corrected antiproton spectra 
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Figure A. 5: Fit parameters for tlie A spectra 
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Figure A. 6: Fit parameters for the A spectra 
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